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_A Generalization of the Method of Conjugate 
Gradients for Solving Systems of Linear 
Algebraic Equations 


It is the purpose of this note to establish a slight generalization of the 
method of conjugate gradients for solving linear equations. The generaliza- 
tion will then be used to provide a unified theory for two rather different 
algorithms which have been proposed for the case in which the matrix of the 
equations is not symmetric. The corresponding formulas for the determinant 
of the matrix will also be developed, although they are probably without 
much practical interest. 

We shall use capital letters for matrices, Greek letters for vectors, and 
lower case Roman letters for scalars. All matrices will be real and of order 
n Xn, and all vectors will be real and of dimensionality . We denote the 
transpose of a matrix or a vector by a prime. 

Let S be a positive definite matrix. The method! of conjugate gradients is 
based upon the construction of two sequences of vectors yo, 71,°--and do, 1, 
--+, and a sequence of scalars do,@;,- --. The sequences are constructed by 
repeated substitution into the following formulas: 


(1) 


= yn — GySén, 


- Y'N4+1YN41 
yN''YN 


(5) by+1 = Yv41 + bvév, N =0,1,---, — 2. 


To avoid discussing exceptional cases, we adopt the convention that if 
by = 0 for some value of N, say N = No, then 6y,41,°--, 5-1 and ywy41,°°°, 
y, are to be defined as being zero vectors. 

HESTENES & STIEFEL' proved a long series of interesting results concern- 
ing the sequence of vectors generated by (1)—(5). In particular, they showed 
that the vectors yo, y1,°°*, Yn—1 are orthogonal (or zero) and the vectors 
50, 51,-- +, dna are S-orthogonal’or conjugate with respect to S. Moreover, 
they show that the scalars ay given by (2) are identical with the coefficients 
in the representation of yo by a linear combination of the vectors Séy. That 


n-l 
Te = 7 a, Sb... 
0 
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From this it can be deduced that the two quadratic forms yy’S-'yy and 
(S-yw)’ (S-*yw) are reduced in value at each step. 

Hestenes & Stiefel also mention without a specific proof that det S = 
(do @1- + *G,_1)~!, where det S stands for the determinant of S. This relation 
is valid only if yo and S are such that dy ~ 0, N < n — 1. A proof of this 
relation can be given as follows: Let A be the matrix whose columns are 4p, 
-++,5,-1, and let G be the matrix whose columns are 7o,° ++, Yn-1. Then by 
the S-orthogonality of the vectors dy, it follows that A’SA = D, where D is 
a principal diagonal matrix whose principal diagonal elements are 40’ So, 
51 S81,° ++, dn-1’Sdn—-1. But from (5) it is seen that éy is the sum of yy and a 
linear combination of the vectors yo, y1,---, yw—1- Therefore by an elemen- 
tary property of determinants, det A = det G. Since the vectors yy are 
orthogonal, G’G is a principal diagonal matrix whose diagonal elements are 
o'Yo, ¥1'71,°**» Y'n-1 Yn—1- Putting these facts together, we get 





det D det D det D 
(7) det S = cua’ detA ~ deta’h ~ detGG 
n—1 §,' S58; 1 





, — . 
k=o Yk Yk Qo 21° * *An-1 


We now consider the applicability of the construction given by (1)-(5) 
to the problem Aé = », where A is an arbitrary non-singular matrix, 7 is a 
given vector, and é is an unknown vector. 

Choose an initial estimate &) of the solution A~y of AE = 7. Choose any 
positive definite S. Let B be an arbitrary non-singular matrix. Construct the 
sequences given by the algorithm with yo = 59 = B(n — Aéo). Then by (6) 


n—1 
B(n - Aéo) = > apSbx, 
0 
or 
n—1 
An = to + © aA“ B-S8,. 
0 
To apply this result, it is advantageous to express S in the form 
S = BATA'B’, 
where T is a positive definite matrix. Then the solution of A-'y is given by 
n—1 
A-'y = f) + > a: TA’B’S,. 
k=) 


The solution can be expressed iteratively, if desired, by the recursion 
relation 


(9) tw41 = Ey + ayTA'B’éy, N= 0,- Ce a= 1. 


From (7), we see that the absolute value of the determinant of A is given 
by 
1 


te © ea an T)4|det B|’ 
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GENERALIZATION OF METHOD OF CONJUGATE GRADIENTS 191 
For some purposes it might be convenient to rewrite the relevant parts 
of the algorithms (1)—(6) in terms of the residual vectors fy = 9 — Aéy. 
To do this we notice that (9) and (3) imply that 
N-1 
B(n ~_ Aéy) = Bon = Boo - > a,BA TA'B’'S; = YN- 
k= 


Thus we have only to replace yy by Boy in (1)—(5) to obtain the desired 
modification. The new relations are 


(1)’ 50 = Bos, 

’ = (Bpn)’ Bon ” 
(2) Ov = 5BATA'B' N= Orem — 1, 
(3)’ Pn+1 = PN — ayATA'B'éy, N= 0,- oo i, 

(Bow+1)' Bows: 

4 i by = ’ N = a" **%, ~~ fy 
(4) IN (Boy) Bow 0 n 2 
(S)’ dv+1 = Bowii + bnin, N =0,---,m — 2. 


These formulas, together with (9), represent the generalization of the con- 
jugate gradient algorithm for the problem A£ = 9 which was promised in the 
first paragraph. 

We shall now specialize T and B so as to obtain algorithms useful in 
practice. 

Case 1. B = I, T = A“, A positive definite. The algorithm (1)’—(5)’, 
(9), assumes the standard fundamental form given by Hestenes & Stiefel." 

Case 2. B = A’, T = (A‘A)—. The formulas (1)’—(5)’ become: 





bo = A'po, 
(A'pw)'A'py 
= oy N=0,---,8—1, 
“N * (Apw)'(Apw) " 
Pv+i = pw — GynAiy, N =0,---,n —1, 
(A’pw+1)’Apn+1 . 
= = — 2 
bw = CA’ ow)’ Aw K-48 


bv41 = A’pyi1 + bwéiy, N =0,---,n — 2. 
The iterations (9) for the solution vector § = Ay are given by 


tyva1 = ty + aviv, N =0,---,n —1, 


and 
|det A| = (do ay: + +d,—1)~?. 


This is the method of solution for the case of an arbitrary non-singular A 
proposed by Hestenes & Stiefel.1 The algorithm constructs a set of A’A- 
orthogonal vectors dy. 
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Case 3. B = I, T = I. The formulas (1)’—(5)’ or (1)—(5) (there is no 
difference in this case) become : 


50 = po, 

Pn'Pn 
y= ON = 0,-+-,n — 1, 
On = §yAA bn’ os 

pn+1 = Pw — GvAA'by, N =0,---,n — 1, 

, 
by = SHON N =0,---,n — 2, 

Pn PN 


ov41 = Pwii + byiy, N =0,---,n — 2. 
The iterations for the solution vector § = A-'y are given by 
fva1 = ty + anA’iy, N =0,---,n —1. 
The formula for the determinant is again 
[det A| = (ao a,- - -a__,)7?. 


The algorithm constructs a set of AA’-orthogonal vectors dy instead of a 
set of A’A-orthogonal vectors, as in Case 2. 

This method of solution for the case of an arbitrary non-singular A was 
proposed by CraiG.? The author learned of this method through a communi- 
cation from Dr. Craig dated May 22, 1953, in which an algorithm essentially 
equivalent to the one given above was stated, but with py calculated at each 
step from the relation py = » — Agy. It has an obvious advantage over the 
Hestenes-Stiefel proposal in that the number of arithmetic operations re- 
quired is slightly smaller. 

To be specific, if we measure the amount of computing work only by the 
required number of multiplications, assuming no zero elements, and count- 
ing a division as two multiplications, the count for the first step of the Craig 
method is 3n? + 5n + 4, whereas that for the Hestenes-Stiefel method is 
4n? + 6n + 4. The detailed breakdown is as follows. First, for the Craig 
method, po is n? multiplications, A’do is then m?, ap is 2n + 2, £; is n, p: is n’, 
p1'p: is n, by is 2 (we already have computed po’po for ao), 5; is ; total 3m? + 
5n + 4. With the Hestenes-Steifel method, po is 2?, 59 is n?, Ado is n?, ao is 
2n + 2, 1: is n, pi is m (it is true that we already have computed Ado, but not 
aoApo nor Adodo), A’p; is n”, (A’p:)'A'p1 is n, 6; is 2, 6, ism; total 4n? + 6n + 4. 
We are obviously assuming full storage of all vectors computed in the previ- 
ous step at each point in the computation. The detailed breakdown is given 
here because the author has encountered a number of apparent discrepancies 
in similar counts published elsewhere (not necessarily relating to the methods 
discussed here); these discrepancies doubtless arise largely from different 
orderings of operations and different assumptions as to what will be stored 
after computation. 

After the first step there is little to choose between the two methods from 
the standpoint of the number of multiplications theoretically required. In 
the Craig method, each step after the first and before the last takes 2m? + 
4n + 4. 


| 





ON 


Ir 
In th 
recur 
it lost 
to ca! 

VW 
for 
numt 
justif 
most 
O(n?) 
error 
speci: 
writi: 
“best 

T 
the a 
be p 
mult 
of A. 
give 
woul 

I 
the r 
(S)’ 
less 1 
the | 
non-! 
in th 
gene! 
does 


New 
New 
1 
syste! 
2 
Mass 
3 


ae 4 


C 


for t 
are | 
met! 
the | 
com; 


bi 
Vy 


mer nN Ds + Oo 


= 


Tt 





ON ESTIMATION OF QUADRATURE ERRORS FOR ANALYTIC FUNCTIONS 193 


In the Hestenes-Stiefel method, the corresponding count is 2m? + 5n + 4. 
In the latter method, if p, in the last step were to be computed from the 
recursion relation, the method would win back the n? multiplications which 
it lost to the Craig method in the first step ; but it would hardly be reasonable 
to calculate the last residual in this way. 

We note that since the algorithm in either case will very probably last 
for n steps unless a peculiarly fortunate choice of & was made, the total 
number of multiplications to be expected is of the order of 2’. The eventual 
justification for the use of conjugate gradient methods in preference to the 
most economical direct methods (which theoretically require only n*/3 + 
O(n?) multiplications) must rest on definitive a priori studies of round-off 
error, and on the peculiarities of individual computing machines, and on 
special characteristics of the matrix A. But it seems quite probable at this 
writing that no final recommendations will ever be formulated as to the 
“best’’ ways of solving linear equations. 

There is a further practical use for the algorithm (1)’—(5)’ other than 
the applications represented by the specializations described above. lt may 
be possible under some circumstances to find a matrix B such that pre- 
multiplication of the equation At = 7 by B will improve the “‘condition’”* 
of A. Then (1)’—(5)’, with T chosen conveniently, say T = (A’A)~, might 
give an advantageous algorithm to use. (With T = (A’A)~ the algorithm 
would be the same as that of Case 2 but with A replaced by B’.) 

In conclusion, we note that Hestenes & Stiefel! assert (p. 424-425) that 
the most general cg-algorithm for any linear system is obtained from (1)’— 
(5)’ by assuming that BA is positive definite and T = (A’B’)-. Neverthe- 
less without this specialization, the generality of (1)’—(5)’ seems to exceed 
the limits specified by Hestenes & Stiefel, because here B is an arbitrary 
non-singular matrix. The disagreement, if it may be called that, seems to lie 
in the fact that their characterization of generality is derived from a certain 
generalization of the method of conjugate directions, and this generalization 
does not include all the possibilities represented by the algorithm given here. 


J. H. Curtiss 
New York Univ. 
New York, N. Y. 


1M. UR. HEsTENES & E. StTIEFEL, ‘Methods of conjugate gradients for solving linear 
systems,’ " NBS Jn. of Research, v. 49, 1952, p. 409-436. 
2E. J. CraiG, Iteration Procedures for ‘Si us Equations. Doctoral dissertation, 
Mass. Inst. Tech., Cambridge, Mass., 1954. 
3 J. Topp, “The condition of certain matrices, I,” Quart. Jn. of Mech. and Appl. Math., 
v. 2, 1949, p. 469-472. 





On the Estimation of Quadrature Errors for 
Analytic Functions 


1. Introduction. The present paper sets forth a complex variable method 
for the estimation of errors which arise when approximate rules of quadrature 
are applied to analytic functions. In contrast to the usual real variable 
methods, this method does not involve the use of the higher derivatives of 
the function, but uses only a knowledge of the size of the function in the 
complex plane. It is therefore of practical value when dealing with situations 
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in which the higher derivatives of the function are difficult to obtain in 
closed form, or where the higher derivatives become excessively large over 
the interval in question. 

The principal feature of this method is the use of the characteristic in- 
equality for bounded linear functionals. The method was developed sub- 
stantially in Davis? and outlined there not only for quadratures but also for 
rules of differentiation, interpolation, etc. Numerical coefficients were ob- 
tained for several rules for the case of functions analytic in the unit circle. 
In the present paper, we consider only quadrature formulas. We deal with 
certain commonly applied rules of low order and also some Gaussian rules 
of rather high order, and present a table of coefficients (Table 2) which per- 
mits the rapid estimation of errors (through formulas (26), (31), (38)) when 
these rules are applied to functions which are analytic on a closed interval 
[a,b]. The estimates given here provide upper bounds for the error and their 
use leads to a conservative practice. In a certain sense, as will be explained in 
§ 4, these estimates are best possible. The problem of obtaining very precise 
estimates of error is a difficult one and such estimates can be obtained only 
after a painstaking analysis of the particular function integrated. When 
estimates are developed which are applicable to wide classes of functions, 
some precision is necessarily lost. It is felt that the present method provides 
a compromise between the time involved in making the estimate and the 
precision of the estimate. It should be added that in any case, the precision 
of the present estimates is independent of the precision of conventional esti- 
mates obtained by using the maximum value of some high derivative over 
the integration interval ; that is, for some functions the former will provide 
the sharper estimate while for other functions the latter will. 

The consideration of classes of functions which are analytic on a closed 
interval (as opposed to functions which are analytic in a fixed circle) pro- 
vides more freedom with but slightly more work. This freedom shows up 
particularly when dealing with functions which are analytic on the segment 
of integration but which have singularities within the circle constructed on 
this segment as a diameter. Functions which are analytic on the interval of 
integration are most conveniently handled by the introduction of the Cheby- 
shev polynomials. These polynomials are orthogonal (with respect to a 
weight) not only over the real interval but over certain ellipses in the com- 
plex plane, and the resulting complex analytic Fourier analysis lends con- 
siderable simplicity to the developments which follow. 

We reiterate that we are dealing with functions which are analytic on a 
closed line segment and hence are continuable analytically to some two di- 
mensional region of the complex plane which contains that segment in its 
interior. Thus, we rule out the cases where the function to be integrated has 
only low order continuity properties or the case where the function has been 
extracted from an ensemble of functions with only certain statistical facts 
known. 

2. The Class L?(&,). Let B designate a bounded region lying in the plane 
of the complex variable z = x + iy. By L?(B), we shall mean the class of 
functions which are single-valued regular analytic in B and are such that 


(1) ff lf (sz) |%dxdy < o, 
B 
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The class L?(B) has been studied extensively. For the fundamentals of this 
theory see BERGMAN.' We note in particular that if f(z) is regular in the 
closure of B, then (1) holds. The positive square root of the quantity (1) 
is generally termed the norm of f over L?(B) and is designated by ||f|| or 
\lfllz. That is, 


(2) ils = (ff eiracay)'s ser). 


B 
If a region G satisfies G C B and if feL?(B), it is clear that feL*(G) and that 


(3) IIflle< Ilfills. 


The principal features of the class L*(B) are as follows. We introduce 
the integral 


Jf see@acay = 0); f, 20°08), 
# 
as an inner product. The class L?(B) possesses complete orthonormal sys- 


tems {f¢,(z)} in the sense of (4). Every function of the class can be expanded 
in a Fourier series 


(5) f(z) = ¥ ang, (2), 
n=0 
wherein 
(6) an = (f, fn) (n = 0,1,---). 


The convergence of (5) is uniform and absolute in every closed subregion of 
B. Moreover, we have 


(7) fl? = ¥ loal?. 


n=0 


The bilinear series }- ¢,(z){,(w) converges for z, weB toa function Kpg(z,), 


a=0 

generally designated as the Bergman kernel function of B, which posseses 
the characteristic reproducing property, (f(z), Ks(z,0)) = f(w) for all 
feL?(B). For simply connected regions B with a Jordan boundary, there 
exists a complete orthonormal set of polynomials {),(z)}. 

There are, moreover, two notable cases wherein the polynomials p,(z) 
have a particularly simple structure. They are the circle and the ellipse. 
The circle |z| < R possesses the set of polynomials 


(8) Pals) = (1 + 1)'e-te*R-"—s (m = 0,1,---) 


as a complete orthonormal system. 

To discuss the case of the ellipse, it is best to assume that it has been 
placed in a normalized position. Let an ellipse, therefore, have its major axis 
along the x axis and its foci at the points (—1,0) and (1,0). Let a and b = 
(a? — 1)# designate its semimajor and semiminor axes respectively and let 
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the quantity p = (a) be defined by 
(9) p= (a+b), a = 3(p§+p-4), b = 3(p* — p-4). 


This ellipse will be designated by &,. For values of p > 1, these ellipses form 
a confocal family and collapse to the segment [—1, +1] asp — 1. In Table 1 
will be found the values of the geometric quantities b and p for a number of 
selected values of a. 

We introduce the Chebyshev polynomials of the second kind by means of 
the definition 


(10) U,(z) = (1 — 2)-'sin ((a +1) arccosz) (nm = 0,1,---). 
It can then be shown that the polynomials 
(11) fa(z) = 2(m + 1)4a-¥(p™*! — p-*"")-1U,, (2), (mn = 0,1,---) 


form a complete orthonormal system for L?(6&,). 

3. Quadrature Errors for L*(&,). In conformity with the above normal- 
ization, we shall assume that the integration to be performed is over the 
interval [ —1, 1 ]. The case of an arbitrary interval may be handled by means 
of an appropriate linear transformation [see § 6]. An arbitrary N + 1 point 
quadrature formula is given by 


+1 N 
(12) [sea ~ E asf) = RU) 


where A, are certain abscissas lying in [—1,1] and a; are the associated 
weights. The error E = E(f) involved in the rule R is 


(13) Ef) = f f(@dx — RU) 


and may be regarded as a linear functional on f. 


If f(x) is analytic on [—1,1] then it is clear that for some value of p > 1, 
f may be continued analytically so as to be regular in the closed ellipse &,. 
Such an f(z) is therefore of class L*(&,), and may be expanded in a series of 
Chebyshev polynomials, 


(14) f(@ = . dat'n(2) : X Jal? = If? < @, 


which converges uniformly and absolutely in the interior of &,. Applying 
the operator E to (14) we obtain 


(15) E(f) = LD a,E ($n). 
n=O 
Estimating E(f) by means of the Schwarz inequality, there is obtained 


(16) IE(N|?< X lanl? E |E(Gs) 
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Let us now write 
ao 
(17) o =o = > |E(f,)|*; 
n=0 


then in view of (7) there is obtained 


(18) IE(f)| < ozllfil. 


The quantity cz which is the norm over L?(&,) of the bounded linear func- 
tional E has an alternate representation in the form o,* = E,E5Ka(z,2). 
It therefore depends only upon the ellipse &,, and the quadrature rule R, 
but is independent of f and of the particular complete orthonormal set of 
functions used in the expansion (14). The quantity og may therefore be 
computed once for all. Using (17) and (11) we have 


(19) at E (n +1) |E(U,(z)) |? 


oa ptt a | ial 
We have, moreover, 
(20) U,(2) = (n + 1)7°T'ngi(z) (nm = 0,1,---), 


where 7,(z) designates the Chebyshev polynomials of the first kind de- 
fined by 


(21) T,(z) = cos (narccosz) (nm = 0,1,---). 


Therefore, 


(22) §:. U,(x)dx = — CT n41(1) 








1 
— Toi(—1)] = GU + (-1) 
If quantities 7, are defined by 
(23) Tm =0 n odd, 
— 
ery n even, 

then we have 

4 C) ( 1 N 2 
(24) oR. = > - a2 ee (+. —_ ; a,U(.)) . 

T n=o P , p k=0 


Table 2 lists the values of « corresponding to the trapezoidal, Simpson’s, 
Weddle’s and the Gaussian 2,3,7,10 and 16-point formulas [For the abscissas 
and weights of these Gaussian rules see Lowan‘ ], and for a range of values of 
the parameter p. These values were computed from (24) on the National 
Bureau of Standards Eastern Automatic Computer (SEAC). 

4. On the Error Estimate (18). The inequality (18) allows us to estimate 
quadrature errors from a table of values of « and from an estimate of the 
norm ||/||. Let us assume that we are dealing with a fixed rule R with error E. 
The value of the parameter p is at our disposal to a certain extent. For func- 
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tions f(x) which are analytic on [—1,1] there will be a range of values of p, 
1 <p < p*< & for which f(z)eLZ*(6&,). For each such p, the inequality (18) 
is valid. To exhibit the dependence of (18) upon p, we should write 


(18’) IE(f)| < or(o)llflle,, 1 <p < 0%. 


Now |lf|l., = 0 when p = 1 and increases as p increases. On the other hand 
or(p) decreases as p increases. Hence, the best estimate will occur for some 


intermediate value of p. The inequality (18’) can therefore be improved by 
writing 


(25) |E(f)| < min or(o)|\flle, 
1<p< p* 


Let the rule R as well as the ellipse &, be fixed. Then the inequality (18’) 
(which is the characteristic inequality for bounded linear functionals) is a 
best possible one in the following sense: there exist functions g of class L?(&,) 
such that 


(26) |E(g)| = or(o)|Iglle,- 


In general, of course, the error will be less. There are two additional facts 
which in the practical carrying out of this method tend to diminish the pre- 
cision of the method. The first is that the exact minimum in (25) is difficult 
to ascertain and what one does is simply to take the minimum of a finite 
number of such values. The second is that the norm of f cannot be computed 
exactly except in the simplest cases and an appropriate upper bound must be 
used. The net effect of all this is to replace the right hand side of (25) witha 
less precise upper bound, but one which is much more readily ascertained. 

5. Methods for Estimating ||f||. It appears, then, that the principal task 
confronting the numerical analyst when using the present method is that of 
estimating the norm ||f|| over some ellipse &,. In the present section we shall 
explain a number of devices which may be useful for this purpose. By the 
very definition (2) we have 


Qn thee = ff se lraxdy, 


& 


and it may be possible in certain simple cases to evaluate (27) directly. If 
f(z) is continuous in the closed ellipse &,, and if we set 


(28) M, = max |f(2)|, 


ze&p 
then we have from (27), 


(29) Iles mt f f axds, 


&p 


so that 


(30) Iflle, < Vad M,, 
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where the quantities a and b are related to p by means of (9). The quantity 
M, or an upper bound for it can, in many cases, be obtained by algebraic 
manipulation. It may be more convenient to replace & by some circle con- 
taining it, say the circle C,: 





p+1 
|z|<a= 2p 
We then have 
(31) IIflle, < IWfllee < ave man \f(z)|, 


assuming that f(z) is regular in Cq. 
If f(z) is of class L?(C,) and its Taylor expansion is known: 





(32) f@) = Xi ans, 
n=0 
then it follows from (7) and (8) that we have 
C) 272n 
(33) Ifllte. = xa? 


which is an exact evaluation of the middle term of the inequality (31). The 
quantity ||f||,, may itself be expressed as an infinite series involving the 
Taylor coefficients a,, but such expressions are more cumbersome and con- 
sequently less useful. 

It should be pointed out that in cases of great complexity, there is always 
the possibility of obtaining M, directly by evaluating |f(z)| along the 
boundary of &, on a sufficiently dense set of points. Despite the loss of ac- 
curacy involved in combining (18) and (30), errors obtained in this way will, 
in general, be better than those obtained by using the conventional real 
variable error expressions and estimating the derivatives which occur there 
by means of Cauchy’s inequality. This can be explained by the fact that the 
real variable expressions must be valid for a wider class of functions. 

6. The Case of an Arbitrary Interval. The quantities of Tables 1 and 2 
refer to the interval [—1,1]. The case of an arbitrary interval [a,b] is dealt 
with by means of the linear transformation 








7 2w b+a 
' bese bw 
(32) _b-a b+a 
w= 2 x+ 2 


which carries the interval a < w < 6 into —1< x< 1. Let the rule E* be 
given on [a,b] as follows: 


b N 
(33) Ev) = [fwd - X asf). 


The analogous rule on [—1,1] is given by 


(34) BY) = [sera - (575) ov [(5) - (F8)] 
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TABLE I: Geometric Quantities for Ellipses 


a b p (wab)! 
1.01 1418 1.3266 .6607 
1.02 .2010 1.4908 .8026 
1.03 .2468 1.6302 8936 
1.04 .2857 1.7574 .9661 
1.05 .3202 1.8773 1.0277 
1.10 4583 2.4282 1.2584 
4.25 .5679 2.9511 1.4324 
1.20 .6633 3.4720 1.5814 
1.25 .7500 4.0000 1.7162 
1.30 .8307 4.5397 1.8420 
1.40 .9798 5.6634 2.0759 
1.50 1.1180 6.8541 2.2953 
1.75 1.4361 10.1515 2.8099 
2.00 1.7321 13.9282 3.2989 
2.50 2.2913 22.9565 4.2421 
3.00 2.8284 33.9706 5.1631 
4.00 3.8730 61.9839 6.9763 
5.00 4.8990 97.9898 8.7723 


and is known by the same name. If f(w) is analytic on [a,b], then 


(35) ee) =s|(25*) 24244] 


2rk b+a h 
oe iva ave 











is analytic on [—1,1] and, setting x, = 





b N a 1 
36) Er) = f" swraw — Eason ="F*| [™ gterae 


2 
~ Beato(52,)] 
so that from (36), (35), and (34), 








(37) ep) = (*5*) zo. 
Thus, 
b— b— 
(8) (Et) = 25 121 < (25*) onlelle, 


The cr(p) are the tabulated values [Table 2] in the z = x + iy plane and 
\|gll, also refers to this plane. 


7. Examples. 


1 
1. Estimate the error E incurred in evaluating f exp (e*)dx by 
0 


Weddle’s rule. From (35) we have g(z) = exp (e+?) which is an entire 
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function of z and is therefore of class L?(6,) for all p > 1. Now, |exp (e(*+”) | 
= exp [Re(e“+”)] = exp [e“*”.? cos y/2]. Thus on &, we have |exp 
(e*+/2) | < exp [e+/?]. By (38) and (30) we have 


(39) |E| <4 (mab)? exp [e@+0?2]g,,. 


According to (18’), the estimate (39) is valid for all a > 1. Of the values 
tabulated, the right member of (39) is minimized for a = 2.5 and yields 


(40) |E| < } (4.242) (315.7) (2.836 X 10-*) = .0019. 


In the above work, the norm of exp (e“+”/?) has been estimated crudely, and 
a number of improvements suggest themselves. Thus, for instance, |exp 


(e(*+0/2)| < exp (e+), and since this last function is concave upward, we 
have 


(41) exp (e@+?) < 3 [exp (e¢+/2) + exp (ee+)/2)] + = [exp e(@+n/2) 
— exp (eet) /2) ], —@¢ 7 x< a. 
It is now easily verified that for arbitrary A,B,C, 
A’a? Bp? 
(42) ff (Ax + By + C)*dxdy = wa (A 5 FP c), 


&p 


so that we have 


A’a? i 
(43) || <4 (rat) (“7 +c) on 
where 
1 
(44) A a 2a [exp (eet) /2) = exp (eet /2) ], 


C = 3 [exp (e(tD/2 a exp (e-e+0/2) ], 
For a = 2.5, we have A = 162.82, C = 158.65, so that 
(45) |E| < } (4.2421) (160.6) (2.836 X 10-*) = .00097 


A conventional estimate of error using the formula 


|E| < max 5 (&) + SF fe @|: f(x) = exp (e*) 
=o tte 1400 
yielded the value | E| < .006. 


4 
2. Estimate the error Z incurred in evaluating f T'(w)dw by means 
3 


of the 7-point Gaussian rule. This is a case where conventional methods can- 
not be used due to the lack of information about the higher derivatives of 
the integrand. Transferring to the interval [—1,1] we must consider the 
function g(z) = I'(4(z + 7)). This function is regular in |z| < 7; hence we 
may take a in the range 1 < a < 7. Now,® |I'(x + iy)| < I'(x) forx >0 
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so that |g(z)| = |T'(}(« + 7) + 34y)| < T[4(@ + 7)]. By the concavity of 
the I function, |g(z)| < max {f'(4(@ + 7)), '(4(—a + 7))}, ze8,. Thus we 
have from (38) and (30), 


(46) |E| < }(wab)* max {T(3(a + 7)), P(3(—a@ + 7))} ce, <a). 
The selection a = 5.0 yields 
(47) |E| < $(8.772) (120) (3.867 10-5) = 2.04 x 10-”. 


P. Davis 
P. RABINOWITZ 
NBS 


ao My a, Kernel Functions and Conformal Mapping, Math. Surveys No. 5, New 
ork, 1950. 

?P. Davis, “Errors of numerical approximation for analytic functions,” Jn. of 
Rational Mech. and Analysis, v. 2, 1953, p. 303-313. 

7C. Lanczos, Introduction, Tables of Chebyshev Polynomials. NBS Applied Math., 
Ser. 3. Washington, D. C., 1952. 

‘ALN. LowaNn, Norman Davins & Artur Levinson, “Table of the zeros of the 
Legendre polynomials of order 1-16 and the weight coefficients for Gauss’ mechanical 
quadrature formula,”” Am. Math. Soc., Bull., v. 48, 1942, p. 739-743. 

, ma e.g., W. Macnus & F. OBERHETTINGER, Formeln und Sdtze, 2nd ed., p. 3, ist 
ormula. 


Integrals Occurring in Problems of Molecular 
Structure 


In theoretical work on molecular structure based either upon the valence 
bond method or upon the molecular orbital method, the electronic wave 
functions are usually built up from atomic orbitals (AO’s). The calculation of 
most physical and chemical quantities then reduces to the evaluation of a 
number of integrals involving these orbitals. It is customary to adopt 
SLATER! type AO’s, defined by (cf. RUDENBERG?’) 


(is) = (§*/m)het (3d#) = (257/3m)tyze-*” 

(2s) = (¢5/m)*re-t (3s) = (2f7/Sx)*(r?/3)-* 
(2pc0) = (¢5/m)4ze-t (3pe) = (267/154) beret 
(2pm) = (£5/x)txe-tr (3px) = (2¢7/15x)*xre-t" 
(2p#) = (f5/x)*ye* (3p#) = (257/15) tyre 

(3do) = (57/2m)*(2? — r°/3)e** = (3d5) = (257/3x) 4h (x? — "Je 
(3dr) = (2¢7/3m)txze-* (3d6) = (2¢7/3x)txye*r 


where (x,y,z) are the cartesian coordinates of the electron referred to the 
nucleus as origin with the z-axis directed along the internuclear axis towards 
the other nucleus, r? = x? + y? + 2 and ¢ is a numerical screening param- 
eter. Occasionally complex Slater type AO’s are used; however since these 
are simply linear combinations of the real Slater AO’s, no distinction will be 
made between integrals involving real and complex orbitals. 

Of the distinct types of two-centre integrals that arise, six are of major 
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importance ; if xa, xa’ etc. denote orbitals of nucleus a and xs, x» etc. orbitals 
of nucleus }, these are* 
Overlap integrals: 


(xa| xs) = f Xaxo dv. 


Direct Coulomb attraction integrals: 


(xa]1/rs|xa) = J xolt/r)xe dv. 


Exchange Coulomb attraction integrals: 


(xe]1/ro|x0) = J xolt/radne a. 


Coulomb repulsion integrals: 


(xara xara) = ff xa(1)x0(2)(1/ras)xw (1) xx (2)dordos 


Hybrid or ionic integrals: 


(xeXa’| Xa’"Xb) = ff Xa(1)xa (2) (1/r12)xa (1)x0(2)doidv2. 


Exchange integrals: 
(xox | Xb’Xa’) _ ff xa(1)x»(2) (1/ri2) xe (1) xa’ (2)dv,dv2. 


in which 7 is the distance of the electron from nucleus ), (1) and (2) represent 
all the coordinates of electrons 1 and 2 respectively and 1:2 is the distance be- 
tween these electrons. The internuclear distance R enters all the integrals as 
a parameter. 

There are essentially only two methods of evaluating such integrals, one 
of which consists of using expansions of the reciprocals of the electron dis- 
tances and the other expansions of the orbitals. The first method is based 
upon the pioneer work of SuGiurA,? ZENER & GUILLEMIN,’ BARTLETT,® 
ROSEN,*® JAMEs’ and others whilst the second is largely due to BARNETT & 
Cou tson® following a suggestion by CooLinGE.? The choice of method for any 
particular problem is somewhat arbitrary and will probably be decided in the 
future by the availability of suitable tables. 

The first method as used in the past has necessitated the computation of 
a large number of auxiliary integrals, defined as follows: 


x +1 
A, (x) -f ein dt; B,(x) -f e-=¢ 2” dt; 


—1 
+1 
G,’(n;x) = f e-** P,?(t)em(1 — £)"de, 
—1 
where P,’(t) is the associated Legendre function of the first kind”; 


felix) = [~ Qe(der**e a, 





wh 


tals 
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where Q,’(¢) is the associated Legendre function of the second kind ;* 


~o +1 
Fr, m,n(%,¥) = f dt f dt! (P — 1)*(E + t)-Oemtine ete ; 
1 —1 
ca +1 
Jn(x,y) = f dt f _ HE MY — 2) — 2) + treme; 
« +1 
K,(x,y) = f dt f dt! (@ — 1)(1 — £2)(1 + w(t + ete"; 
1 -1 
«a +1 
Frm, n(%,9) -f dt f de! (t + ¢)-*(— #9(1 — 8) ete"; 
1 —1 


C) t 
S,(m,n;x) = f dtQ,(t) e~** a f e~=*' g's dt’; 
1 1 
o;(m,n;x) = f dtQ, (t) e~** a f est gs dt’: 
1 t 
H, (m,n ;x) = S,(m,n;x) + S,(n,m;x); 
W'(m,n;x) = f dt f dt'Q,?(t >)Per(t <)eree+egmy'n (g'2 — 1)ri2, 
1 1 


All of these auxiliary integrals are expressable in terms of polynomials, ex- 
ponentials, natural logarithms and the exponential integral —Ei(— x) = 


f te-‘dt. Many of them are simply related to each other. There is a 


x 
distressing lack of uniformity not only in the choice of auxiliary integrals 
computed but also in the notation. 

The second method, which is particularly valuable in the case of integrals 
involving more than two centres, requires the computation of certain so- 
called ¢-functions which are expressable ultimately in terms of Bessel func- 
tions of half-integral orders for which the argument is pure imaginary." 

Several approximations to the exchange and hybrid two-centre integrals 
and to multi-centred integrals have been suggested” and their accuracy has 
recently been assessed by comparison with the exact values". 

This review concludes with three bibliographies ; 


A—Reduction of two-centre molecular integrals 
B—tTables of two-centre molecular integrals 
C—Miulti-centre integrals. 


A very large number of papers relevant to molecular integrals has been pub- 
lished and it is not possible (nor desirable) to describe all of them. Bibliog- 
raphy A is a somewhat arbitrary selection of papers chosen so that any 
integral that has been calculated to date may be found in at least one of 
them and chosen to include also those papers which the reviewer considers 
to be the most useful. Bibliography B describes in as much detail as space 
allows the tables of two-centre integrals and of auxiliary functions which 
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have been published since the date of inception of MTAC in 1943. Excluded 
from this bibliography is the multitude of published material which contains 
only a few isolated integrals and whose usefulness is effectively restricted to 
the problem concerned in the particular reference. Warning must be given 
that owing to the complexity of the work, the occurrence of errors is frequent, 
even the corrections being by no means exhaustive or completely reliable. 
The final bibliography C lists without description (apart from title) those 
papers concerned with multi-centre integrals. No systematic tabulation of 
multi-centre integrals has yet been published. 

Information on tables published prior to 1943 may be obtained from the 
FMR Index" and from the excellent bibliographies of MuULLIGAN,” Root- 
HAAN® and RUDENBERG."” 

Finally it seems worthwhile to mention that a discussion of future de- 
velopments was reported in the proceedings of a conference on quantal 
methods in valence theory (QMVT).!8 


Bibliography A 
Reduction of two-centre molecular iniegrals. 


N. Rosen, “Calculations of interactions between atoms with s-electrons,”’ 
Phys. Rev., v. 38, 1931, p. 255-276. Almost all types of molecular integrals 
involving identical nus orbitals. 


M. Kotani, A. AMEmMiya & T. Simose, ‘Tables of integrals useful for the 
calculations of molecular energies I,’’ Phys.-Math. Soc. Japan, Proc., v. 20, 
1938, Extra no. 1, p. 1-70. Most of the integrals involving = 2 orbitals 
with a common screening parameter and a large number over 1s and n = 2 
orbitals, the » = 2 orbitals having a common screening parameter. 


C. A. Coutson, ‘“Two-centre integrals occurring in the theory of molecular 
structure,” Camb. Phil. Soc., Proc., v. 38, 1941, p. 210-223. A wide selection 
of one electron integrals over ms and npo orbitals with no restriction on the 
screening parameters. 


R. G. Parr & B. L. CRAWForD, JR., “On certain integrals useful in molecular 
orbital calculations,” Jn. Chem. Phys., v. 16, 1948, p. 1049-1056. Most 
types of integrals over like 2p7 orbitals. 


R. S. MuLLIKEN, C. A. RIEKE, D. Ortorr & H. Ortorr, ‘Formulas and 
numerical tables for overlap integrals,” Jn. Chem. Phys., v. 17, 1949, p. 
1249-1267. Overlap integrals over all combinations of ns, npo, and npx 
AO’s for » = 1,2,3,5 with no restriction on the screening parameters. 


H. J. Koprneck, “‘Austausch und andere Zweizentrenintegrale mit 2s- und 
2p-functionen,”’ Zeits. Naturforschung, v. 5a, 1950, p. 420-431; ‘‘Zweizen- 
trenintegrale mit 2s- und 2p-functionen II. Ionenintegrale,’’ Zeits. Natur- 
forschung, v. 6a, 1951, p. 177-183. All the integrals involving only m = 2 
orbitals given in the second reference and some additional hybrid integrals. 


G. Arakr & W. Watart, “Electronic states of C2-molecule I1[I]—Numerical 
values and reduction formulae of integrals,” Prog. Theoretical Phys., v. 6, 
1951, p. 961-979. The hybrid integrals over like m = 2 orbitals not given in 
reference (A2). 
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R. O. BRENNAN & J. F. MULLIGAN, “‘Two-center heteronuclear hybrid 
coulomb-exchange integrals,” Jn. Chem. Phys., v. 20, 1952, p. 1635-1644. 
All hybrid integrals involving = 2 orbitals with the restriction that orbi- 
tals on the same atom have a common screening parameter. 


C. C. J. RoorHaan, “A study of two-center integrals useful in calculations 
on molecular structure I,” Jn. Chem. Phys., v. 19, 1951, p. 1445-1458. 
All molecular integrals other than hybrid and exchange integrals over all 
combinations of orbitals with 2 = 1 and 2 with no restriction on the screen- 
ing parameters. 


K. RODENBERG, “‘A study of two-center integrals useful in calculations on 
molecular structure II. Exchange integrals,” Jn. Chem. Phys., v. 19, 1951, 
p. 1459-1477. Exchange integrals over all combinations of orbitals with 
n = 1,2,3 with no restriction on the screening parameters. 


M. Korant, E. IsnHicuro, K. Hijrkata, T. NAKAMURA & A. AMEMIYA, 
“Tables of integrals useful for the evaluation of molecular energies III,” 
Phys. Soc. Japan, Jn., v. 8, 1953, p. 463-476. All the integrals over like 
n = 2 orbitals not given in the earlier papers in this series. 


These references are based upon the first method. The remaining refer- 
ences use the second method. 


M. P. Barnett & C. A. Coutson, “The evaluation of integrals occurring 
in the theory of molecular structure,” Roy. Soc. Phil. Trams., v. 243A, 
1951, p. 221-249; “The evaluation of unit molecular integrals,” QMVT, 
1951, p. 238-271. All molecular integrals. 


Bibliography B 
Tables of two-centre molecular integrals. 


V. GriFFINnG, “The size and vibration frequency of the excited benzene 
molecule,” Jn. Chem. Phys., v. 15, 1947, p. 421-430. Most types of integrals 
over like 2pz orbitals for 2R = 8-87, 8-67 (0-1) 8-17 (0-2) 7-77. 


R. G. Parr & B. L. Crawrorp, Jr., “Molecular orbital calculations of 
vibrational free constants I. Ethylene,” Jn. Chem. Phys., v. 16, 1948, p. 
526-532. Many integrals over like 2px orbitals for 25R = 6, 8-13, 8-37, 10. 
“On certain integrals useful in molecular orbital calculations,’’ Jn. Chem. 
Phys., v. 16, 1948, p. 1049-1056. Many integrals over 2px orbitals for mainly 
2¢R = 0, 3 (1) 10. 

R. S. MutiiKen, C. A. Rreke, D. Ortorr & H. Ortorr, “Formulas and 
numerical tables for overlap integrals,’ Jn. Chem. Phys., v. 17, 1949, p. 
1249-1267. Overlap integrals over all combinations of m = 1, 2, 3, 5 orbitals 
usually for (¢ + ¢’)R = 0 (1) 16 (2) 20, (§ —*)(¢+%") = —-6 (-1) -6. 
Also some linear combinations of these. 

P. J. WHEATLEY & J. W. Lrnnetrt, “Molecular force fields Part XI—A wave 
mechanical treatment of the change with distortion of the interaction energy 
of carbon 2px orbitals,’ Farad. Soc., Trams., v. 45, 1949, p. 897-902. Most 
types of integrals over like 2px orbitals for [R = 3-5_(0-5) 5. 
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J. O. HigscHFELDER & J. W. LinneETT, “The energy of interaction between 
two hydrogen atoms,” J. Chem. Phys., v. 18, 1950, p. 130-142. So(m,n,x) for 
x = -5 (-5) 5 (1) 8(2) 12, 1-75; m = 1 (2) 5, n = 0 (2) 6. H(m,n,x) forx = 
0-5 (0-5) 5 (1) 8(2) 12, 1-75; m = 0, 2, nm = 0 (2) 6 and m = 4, n = 4, 6 
and m = 6 = n. Many molecular integrals over 1s and ” = 2 orbitals, all 
orbitals having same screening parameter for ¢R = 0, 1-5, 1-75, 2 (-5) 5 
(1) 8 (2) 12. 


H. J. Koprineck, ‘‘Austausch- und andere Zweizentrenintegrale mit 2s- 
and 2p-functionen,” Zeit. Naturforschung, v. 5a, 1950, p. 420-431. Most 
molecular integrals except hybrid integrals over » = 2 orbitals with a com- 
mon screening parameter ¢ for {R = -5 (-5) 7 with a few omissions. 


“Zweizentrenintegrale mit 2s- und 2-Funktionen II,” Zeit. f. Naturfor-. 
schung, v. 6a, 1951, p. 177-183. Hybrid integrals over m = 2 orbitals with a 
common screening parameter ¢ for ¢R = -5 (-5) 7. 


“‘Zweizentrenwechselwirkungsintegrale III. Integrale mit 2p- und wasser- 
stoffanhnlichen 2s-Funktionen,” Zeit. Naturforschung, v. 7a, 1952, p. 785- 
800. All types of integrals over 2p and a linear combination of 1s and 2s 
orbitals with common screening parameter ¢ for ¢R = 1 (-5) 7. 


G. Araki & W. Watar!I, “Electronic states of C2-molecule,”’ Prog. Theor. 
Phys., v. 6, 1951, p. 961-979. Many integrals over m = 2 orbitals with a 
common screening parameter ¢ for ¢R = 1 (-25) 5. 


M. P. Barnett & C. A. Coutson, “The evaluation of integrals occurring in 
the theory of molecular structure,” Roy. Soc. Phil. Trans., v. 243A, 1951, 
p. 221-249. Bessel functions of the first and second kinds of pure imaginary 
argument In43(x) and Kn44(x) for m =— 1 (1) 4,x = 0 (-5) 10 and K,4;(x) 
for x = 10 (1) 25. 
A. A. Frost & J. BRAUNSTEIN, “Hydrogen molecule energy calculation by 
correlated molecular orbitals,”” Jn. Chem. Phys., v. 19, 1951, p. 1133-1138. 
Fo10(2x, x) =i Foo1(2x, x), Foxo(x, 0) = Fooi(x, 0) and Fo10(x, _ x) — Foo 
(x, —x) for x = 1 (-5) 5 (1) 7 and certain special integrals over like 1s 
orbitals for ¢R = 0 (-5) 5 (1) 7. 
K. F. HERZFELD & V. GriFFING, “Work on integrals at the Catholic Uni- 
versity of America since 1949,”” OMVT, 1951, p. 202-204. Direct Coulomb 
attraction integrals over 2p, 3p and 4p orbitals with a common screening 
parameter for 2¢R = 5, 8, 10 (5) 20. 
E. IsHicuro, K. Hijikata, T. Arat & M. Mizusnima, Natural Science 
Reports, Ochanomizu University, v. 1, 1951, p. 22. Certain types of auxiliary 
integrals (cf. M. Korant,!* 1951, p. 210). 

[The reviewer has been unable to obtain this paper. ] 


J. F. Mutuican, ““LCAO self-consistent field calculation of the ground state 
of carbon dioxide,”’ Jn. Chem. Phys., v. 19, 1951, p. 347-362. Many integrals 
over 1s(¢) and m = 2(¢’) orbitals for ¢ = 5-70, 7-70, ¢’ = 1-625, 2-275 and 
R = 2-2. 

R. S. MuLLIKEN, “Overlap and bonding power of 2s, 2p-hybrid orbitals,” 
Jn. Chem. Phys., v. 19, 1951, p. 900-912. Overlap integrals for 2s, 2p-hybrid 
AO’s of like atoms for ¢R = 3-0 (0-2) 5 (0-5) 7. Some additional overlap 
integrals. 














———— 
-~ ——— 
NS 












t 2s- 
Most 
com- 


rfor-. 
ith a 


sser- 
785- 
d 2s 


heor. 
th a 


ig in 
951, 
lary 
4(x) 


1 by 
138. 
Foo 
e Is 


Tni- 
ymb 
ling 


nce 
ary 


‘late 
rals 
and 
ls,”’ 
rid 
lap 





INTEGRALS OCCURRING IN PROBLEMS OF MOLECULAR STRUCTURE 209 


E. Scrocco & O. Satvettt, “Sul calcolo di certi integrali utili nel metods 
degli orbitale molecolari antisimmetrizzati,” La Ricerca Scientifica, v. 21, 
1951, p. 1629-1639. (2p. 2pxa:|2px. 2pm) for ¢R = 0, 4 (2) 10, 14, 20; 
t’/t = -6 (-2) 1-4 with a few omissions. 

“Studio di alcuni integrali che si presentano nei problemi di struttura 
molecolare,”” La Ricerca Scientifica, v. 22, 1952, p. 1766-1775. (2px. 2pme-| 
2pxa 2pm) for ¢R = 4 (2) 10, 14, 20; ¢’/¢ = -6 (-2) 1-4; 


“Studio dell integrali di scambio e di ricerca di formule approximate per 
alcuni integrali che si incontrano nei problemi di struttura molecolare,” 
La Ricerca Scientifica, v. 22, 1953, p. 98-107. (2p%. 2pmy|2pxy 2px_) for 
tR = 4 (2) 10, 14, 20; ¢’/¢ = 0-6 (0-2) 1-4 


Y. Aoxt, “Numerical tables of overlap integrals for 3s, 3p hybrid A — 0,” 
Phys. Soc. Japan Jn., v. 7, 1952, p. 451-457. Overlap integrals between 1s, 
2s, 2pe(¢) and hybrid (3s, 3pc) orbitals (¢’) for (¢ + ¢’)R = -0 (1) 6 (-4) 
10 (1) 16 (2) 20 and (¢’ — £)/(¢’ + $) =— -4 (-1) 0-3; between hybrid (3s, 
3pc) orbitals for ¢’R = 0 (5S) -2 (8) -5 (10-0). 

R. O. Brennan & J. F. MuLuican, “Two-center heteronuclear hybrid 
coulomb-exchange integrals,”’ Jn. Chem. Phys., v. 20, 1952, p. 1635-1644. 
Hybrid integrals over n = 2 orbitals such that orbitals on same atom have a 
common screening parameter for ¢ = 1-625, ¢’ = 2.275, (R = 2.2). 

I. FiscHer, “The approximations of molecular theory tested on LiH and 
BeH+,” Arkiv for Fysik, v. 5, 1952, p. 349-376. Some integrals over 1s 
(¢ = 1, 2-69) and 2s(¢ = -658) and 2p(¢ = -545) orbitals for R = 2-64, 
3-02, 3-40 and 1s(¢ = 1), 2s(¢ = -979), 2p(¢ = -871) for R = 2-48. 


2 
C. W. Jones, A Short Table for the Bessel Functions I,.4(x), = Kata(*), 
1952, University Press, Cambridge. x~*~! J,4;(x) and = ett Kas;(x) for 


0 (-1) 5, 2 = O (1) 10; e-* Tn44(x) and Se Kasi(x) for x = 5 (-1) 10, 
0 (1) 10. 

T. Murat & G. Arak, “‘Calculation of heteronuclear molecular integrals,” 
Prog. Theor. Phys., v. 8, 1952, p. 615-638. J,(x,y) and K,(x,y) for n =— 3 
(1) 2 and x, y usually integers < 10 with some additional half integral values. 
Fr, n(x,y) for m = 1, n = 0 (1) 9, x = 1 (1) 6 (1) 10, y = 0; also a large 
number of values of Fn,.(x,y) for m,n integers and x, y usually integers in 
the ranges O< m< 4, O< 2» < 9, 25 x< 10, 0-S< y< 6-5. Integrals, 
mainly coulomb and hybrid over m = 2 orbitals such that any two orbitals 
on the same atom have a common screening parameter. 


x 


n 


I. G. Ross, ‘‘Calculations of the energy levels of acetylene by the method of 
antisymmetric molecular orbitals, including o — # interaction,” Farad. 
Soc., Trans., v. 48, 1952, p. 973-991. Many integrals over » = 2 orbitals 
with a common ¢ for ¢{R = 3-61. 


M. Koran, E. Isaicuro, K. Hrjyrkata, T. NAKAMURA & A. AMEMIYA, 
“Tables of integrals useful for the calculations of molecular energies III,” 
Phys. Soc. Japan, Jn., v. 8, 1953, p. 463-476. Extensive tables of molecular 
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integrals over = 2 orbitals having a common screening parameter for 
¢R = 2 (-25) 7. 

H. H. Jarré & G. O. Doak, “On the stability of metallo-organic compounds,” 
Jn. Chem. Phys., v. 21, 1953, p. 196-200. Overlap integrals over 2s, 5s, 2c, 
Spo and 3de orbitals mainly for (¢ — ¢’)/(¢ + ¢’) =— -75 (-25) -75, 
(¢ + ¢’) = 0 (1) 16, 18, 20 (4) 32. 


H. H. Jarré, “Some overlap integrals involving XX orbitals,” Jn. Chem. 
Phys., v. 21, 1953, p. 258-263. Overlap integrals over 2p2, 3dx, 5dx orbitals 
for (¢ — £')/(§ + £’) = —-75 (-25) -75, (¢ + ¢’) = 0 (1) 16, 18, 20 (4) 32. 
C. J. THorNeE, R. S. BARKER & H. Eyrine, ‘“‘Tables of quantum mechanical 
integrals. I. Some two parameter integrals,’’ Studies in Applied Math. no. 10 
(Univ. of Utah, Salt Lake City), 1953. Overlap and Coulomb repulsion 
integrals and others over 1s orbitals for R = 1 (-1) 1-5 (-05) 2 (-1) 5;¢ = 
-7 (-05) 1 (-02) 1-5. 

D. P. Craic, A. Maccoty, R. S. Nynoim, L. E. Orcer & L. E. Sutton, 
“Chemical bonds involving d-orbitals. Part II,”” Chem. Soc. Jn., Jan. 1954, 
p. 354-357. Overlap integrals over some of the possible combinations of 
n = 3 and m = 4 orbitals usually for (¢ — ¢’)/(¢ + ¢’) =— -5 (-1) <5, 
(¢ + ¢R = 4 (2) 12 or 20. B,(x) for x = 0 (-1) 1-5, m = 0 (1) 8. 

R. S. BARKER, H. Eyrine, C. J. Toorne & D. A. Baker, “Use of electron 
repulsion integral approximations in molecular quantum mechanics,” Jn. 


Chem. Phys., v. 22, 1954, p. 699-702. Hybrid and exchange integrals over 1s 
orbitals for R = 1 (-05) 3 (1) 6, 8, 12. 
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Characters of the Symmetric Groups of 
Degree 15 and 16 


Introduction. The construction of the complete character table for the 
symmetric group of a given degree n (i.e., the group formed by the permuta- 
tions of m quantities) is a problem of long standing. Theoretically the prob- 
lem was solved many years ago by FROBENIUS,! who first introduced the con- 
cept of the character of a representation (thus meriting the title of founder 
of the theory of group representations). Frobenius’ expression for the sym- 
metric group characters as coefficients of a certain algebraic form is, however, 
unsuitable for practical calculation except for groups of quite low degree. 
Some years ago MURNAGHAN? was able to derive a recursion scheme which 
permits one to calculate the characters for a given m in terms of the (pre- 
sumably known) characters of all the symmetric groups of lower degree. 
This scheme has been employed by various authors* to construct character 
tables for all the symmetric groups of degree » < 14. The recursion scheme 
effectively reduces the problem to an exercise in bookkeeping, but, despite its 
simplicity, it suffers from two practical faults. First of all, the labor involved 
in computing the characters by hand becomes prohibitive for m larger than, 
say, 14. Secondly, as mentioned above, for each m the results for all lower a 
must be used, so that errors may be expected to propagate. This, of course, 
is not a theoretical limitation, since in principle each table may be calculated, 
ab initio, and any given table can be checked for consistency by using the 
fact that the characters themselves satisfy certain orthogonality relations. 
The impracticability of proceeding by hand computation is clear when it is 
remembered that the number of characters of the symmetric group of degree 
m is equal to the square of the number of unrestricted partitions of ; for 
n = 15, this number is 30,976. 

It should be evident that the practical difficulties mentioned in the pre- 
ceding paragraph can in large measure be overcome by the use of electronic 
computers. In fact, it seemed to us that this very calculation would provide 
a very good test of the speed and flexibility of such a computer with regard 
to the handling of purely algebraic problems. Accordingly, we set up and 
carried out the calculation of the complete character tables of the symmetric 
groups of degree m = 15 and 16 using the Los Alamos Electronic Computer 
(MANIAC). The method was checked by re-computing the tables‘ for 
n = 10 tom = 14. The method of calculation is described in Sec. (III). The 
experience was quite encouraging, and suggests that it would be profitable 
to apply electronic computer techniques to a large class of quite complicated 
problems in algebra and group theory. 

The Character Formulae. For convenience we give here a brief resume of 
the relevant properties of the symmetric group characters. For proofs and 
detailed discussion the reader is referred to the standard works of LITTLE- 
woop‘ and MuRNAGHAN.® 
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The basic fact concerning the symmetric group of degree m is that both 
its classes and its irreducible representations are in one-to-one correspondence 
with the number of unrestricted partitions of the integer . More precisely, 
each class may be denoted by a partition of m into non-negative integers: 


(p) = (p1, p2- - *pe) 
(1) pir p2 2--*2D Oe 


k 

Le=n 

i=1 

Similarly, each irreducible representation may be so labelled. In the 

sequel we shall reserve the symbol (p) to label classes, and the symbols (A), 
(u) to label irreducible representations. In the interest of compactness we 
adopt the standard convention of denoting repeated parts by a super-script ; 
thus the partition (4,3,3,3,1,1) of 15 will be written (4,3*,1?). With regard to 
the labelling of classes, this notation has a simple interpretation in terms of 
the one-dimensional permutation representation of the group, the class 
(4,3*,12) being that which consists of one cycle of order 4, three cycles of 
order 3, and two cycles of order 1. For a given irreducible representation 
(X), the character of the class (p) (i.e., the common trace of the matrices be- 
longing to this class and representation) will be denoted by x}. For the 
symmetric group, x is always an integer, positive, negative, or zero. If we 
consider the partition (p) as labelling the rows and (A) as labelling the 
columns, the quantities x?) are the elements of a matrix of degree equal to 
the number of partitions of ». The x2} satisfy the following orthogonality 
relations: 


(2a) , * ex Dx = niin 
(pe) 
n' 
r ! 
(2b) ) xx, = TF Sop 
@) h, 


where h, is the number of group elements belonging to the class () and is 
given explicitly by 
(3) h, 


n! 
~ Jaq! 288! 37y!--- 





a being the number of unary cycles in (p), 8 the number of binary cycles, etc. 
For every class, the character of the irreducible representation (m) (cor- 
responding to the partition of m into one part) is unity; i.e., 


(4) x = 1 for all (9) 


Similarly, for every class, the character corresponding to the partition 
(1") is + 1, depending on whether the class (p) is even (i.e., has an even 
number of cycles of even order) or odd. 


5) x” =+1, (p) even 
x =— 1, (p) odd 


If we denote by (X) the partition of m conjugate’ to (A), then for every 
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class, the characters of the irreducible representations corresponding to (X) 
and (A) are related by: 
(6) xO = xx 

The consequence of this last relation is that for each (p) we have to 
calculate only slightly more than half the characters, since the character of 
the associated representation is either equal to, or just the negative of, the 
character of the given representation. If P(m) denotes the total number of 
unrestricted partitions of m and P,(m) the number of self-conjugate parti- 
tions, then the total number* of quantities x?} which must be calculated is 
just 4P(n) [P(n) + P.(n)]. For n = 15 this number is 15,840. 

The recurrence relations for the characters referred to in the introduction 
may be stated in the following form :5 

Let (p’) be a partition of the integer (7 — m) which differs from the 
partition (p) of by the omission of a cycle of order m. Then 


(7) x} = x + x, 
” 


where the partitions (u) of (” — m) run over the set 
(8) (uw) = (Ar — m, do,-+ de), (Ar, A2 — ™M, Age + Ae)y+ ++, (An, Aes + *Ae — m) 


The partitions (8) do not in general satisfy the conditions (1), since the 
parts will not necessarily be in non-ascending order or even positive. It can 
be shown®:* that the following three rules may be invoked: 


1. If (u) hasall its parts in non-ascending order but the last part is negative, 
then the contribution x“ to the sum (7) is identically zero. 
2. If the ith part of (u) is one less than the ¢ + ist part, wy +1 = piyi, 
then the corresponding contribution to (7) is identically zero. 
3. x™ =— x“ where (u”’) is the partition obtained from (u’) by the re- 
placements 
9) ——— 
Mi+1 —> ph 
The procedure is thus to operate on the set (8) with rule 3 until either a 
partition satisfying (1) is obtained or the conditions of rule 1 or 2 are satis- 
fied. (The + sign in (7) is intended to represent possible sign changes arising 
from rule 3.) In conjunction with these operating rules, (7) effectively re- 
duces the calculation of the characters of a given group to those of groups of 
lower degree. By iteration of this process a complete reduction can be carried 
out for every given (p) until (9) = (0). Of course no useful closed expression 
for this process exists. Formulae do exist®:* which give x?) directly in terms 
of the number of cycles of various orders occurring in (p), but they become 
prohibitively complicated unless (A) has a particularly simple form. There is, 
however, one explicit formula, due to Frobenius, which is of use in a practical 
calculation, namely that which expresses the character of the identity class 
(p) = (1”) in terms of the parts of (A). This may be written in the form® 


n! IL (1; => 1;) 


(10) xX = hikeokl (A = Aa, Ag+ ,An) 
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where /; = \; + k — 4. This gives the trace of the unit matrix in the repre- 
sentation (A) and is hence the dimension of the representation.® 

This completes our resume of the properties of the characters of the 
symmetric group. Many other relations exist, but they are not of interest 
for our purpose. 

Calculation Scheme. The computer first simply prepares a table of the 
partitions of ”. This table is also used to designate the classes. In addition 
the complete character tables for the symmetric groups of degree n < 4 are 
stored in the computer. Their function will be described presently 

The sequence of steps performed by the computer is as follows: a character 
is first decomposed by (7) into a sum of characters all of the same class 
(p’), corresponding to a group of degree n’(< m). Attention is fixed on the 
first component. The integers describing the partition are rearranged, if 
necessary, in non-ascending order, with the aid of (9). Next the partition is 
examined to see if 


(i) its value is zero (cf. rules 1 and 2, page 214). If it isso, we proceed to 

the next component of the decomposition ; if not, we test to see if 

(ii) 2’ < 4. If this condition is satisfied, the computer makes reference 
to the appropriate character table stored in the memory; if this 
condition is not satisfied, the computer finally determines if 

(iii) (p’) = (1*’). If so, the character is computed by (10). If not, a 
further reduction by means of (7) is made and the above procedures 
repeated until the first component is evaluated ; all members of the 
original decomposition are treated in turn. 


Obviously, the most efficient method of utilizing (7) is to choose the reduc- 
tion parameter m as large as possible, and this was always done. 

The instructions, or code, for the problem occupied about one-half of the 
1024 locations of the electrostatic memory ; numerical storage accounted for 
the remaining half; the partitions occupied a small fraction of the auxiliary 
drum memory. The computing time for the group » = 15 was five hours; 
for m = 16, twelve hours were necessary to compute the 27,258 distinct 
characters. 

Obviously, the number of operations necessary to compute a character 
is very heavily dependent upon the structure of the partition (p). Indeed, for 
the class (p) = (2°) six levels of branching were required to obtain the nec- 
essary reduction and the computer took thirty-seven minutes to calculate 
the characters for this class. In comparison only one minute was needed to 
compute the characters for (p) = (16). 

Although electrostatic storage limitations precluded a complete orthogo- 
nality check on the table, (2b) was used for each class separately. 

The two tables of characters have been photographed on microfilm di- 
rectly from the computer print-out. Each frame contains the characters of 
one class. These films are available in the UMT Fite [UMT 195]. 


RoBeErtT L. BIvINs 
N. METROPOLIS 
Pau. R. STEIN 


Los Alamos Scientific Lab. Mark B. WELLS 


Univ. of Calif. 
Los Alamos, New Mexico. 
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1G. FrRopentus, Preuss. Akad., Berlin, SitzBer., 1900, p. 516. A list of Frobenius’ 
relevant pers may be found in ref. 4 or 5. 


MURNAGHAN, Am. Journal of Math., v. 59, 1937, p. 437. See also the paper of 
a 2 & RICHARDSON.? 


3D. E. LittLEwoop & A. R. RicHARDsON, Roy. Soc., Phil. Trans., v. 233A, 1934, p. 99. 
4M. Z1a-up-Din, London Math. Soc. Proc (2), v. 39, 1935, p. 200; ibid., v. 42, 1937, 
p. 340. 


K. Konno, Phys. Math. Soc. Japan, v. 22, 1940, p. 585. 
5D. E. LittLEwoop: The Theory of Group Characters, 2nd edition. Oxford, 1950. 
6 F. J. MuRNAGHAN: The Theory of Group Representations. Baltimore, 1938. 


7 The partition conjugate to (A) is obtained by interchanging rows and columns in the 
Ferrers-Sylvester graph of (A). 
® P,(n) = P(n) — 2P(n — 2) + 2P(m — 8) — 2P(m — 18) +---. This was pointed 
out to us in a private communication from Professor N. J. FINE. 
® The following question is of some practical interest : given , for what partition (A) does 
(9) take its largest value, and how does this value vary with n? The authors have been unable 
to find any discussion of this problem in the literature. 
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1234[A].—M. Lorxin & M. E. Younc, Table of Binomial Coefficients. 
Exact Values. Ballistic Research Laboratories Memo. Report No. 762 
Aberdeen Proving Ground, 1954, 49 p., 21.6 X 27.9 cm. Mimeographed 
from typescript. 

This table is a sequel to RMT 1123. It gives exact values of the coeffi- 
cients 
n! 
ri(n — r)! 


for r< (n + 1)/2 and m = 0(1)100, whereas the previous table rounds 
these values to 20 figures. This new table will be of use in studies involving 
congruence properties and other theoretical properties of binomial coeffi- 
cients. 


D. H. L. 


1235[(C,D,E].—R. A. Hirvonen, “Nutshell tables of mathematical func- 
tions for interpolation with calculating machines,” Bull. Géod. No. 30, 
1953, p. 369-392. 
The tables presented here are one page tables. Interpolation is by means 
of Taylor’s development which the author writes in the nested form 


fa+th) = {((Dt+ C++ BY + Ali+f @ 


so that 
A =hf'(a), B= hf''(a)/2!},... 


The tables therefore give the functions at coarse argument intervals together 
with the coefficients A,B,C,D. The tables are as follows: 


In x for x = 1(.02)1.6; 10D 
e* for x = 0(.01).2: 10D 
sin x for x = 0°(2°)90°; 10D 
arctan x for x = 0(.02)1; 8D. 
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The latter table is in decimal parts of a degree. Three special geodesic tables 
are used to illustrate Taylor interpolation. 
D. &. L. 


1236[F ].—H. Gupta, M.S. Coeema & O. P. Gupta, “On Mdbius Means,” 
Panjab Univ. Research Bulletin No. 42, 1954, 16 p. 
The table occupying all but the first page of this work was computed in 
order to study conjectures of V. Brun and C. L. SIEGEL regarding the second 
numerical integral of the Mébius function u(m), the multiplicative function 


for which 
u(p*) = — [1/a]. 
If we define 


ui(m) = Su(k) and a(n) = E mlb) 


k=1 
then the conjecture of Brun is that “on the average” 
po(m — 1) =— 2n + 12 + Bn 
where 8 = — 18. Siegel considered the value 
B =— 2x*/f(3) =— 16.42119. 
The authors tabulate 
F(n) = To(n) = mpeo(m — 1) + 2m? — 12n 


whose average value should be 8, and the successive means T; defined by 
nT,(n) = Y Tir-a(») 
v= 


for k = 0(1)5, m = 1(1)750. Exact values of T)(”) are given. The other T’s 
are given to 3D. Graphs of T;, T, and 7; are given. For 375 < m < 750, T5(m) 
descends nearly monotonely from —14.3 to —15.0. 

D & &. 


1237[K ].—J. H. CALDWELL, “‘Approximating to the distributions of meas- 

ures of dispersion by a power of x’,”’ Biometrika, v. 40, 1953, p. 336-346. 

It is assumed that all observations are from normal populations. If the 
measure of dispersion is u, then for suitably chosen c and X, cu* has approxi- 
mately the x? distribution with v degrees of freedom. Table 1 of this paper 
gives values of v, A, log c for the range and mean deviation for sample size 
n,n = 2(1)20 and for the first quasi-range with m = 10(1)30. vis given to 2S, 
d and log c to 5S. The first quasi-range is the difference between the largest- 
but-one and the smallest-but-one observations in a sample of m. Table 2 
contains the same constants for the mean of m ranges m = 2(1)5 and = 
2(1)10, where 1 is the size of each of m samples. Tables 3 and 4 give the ap- 
proximate upper 5% and 1% points of the ratio of maximum value to mini- 
mum value in a set of k independent ranges (each range for a sample of size 
n from the same population). Values are given to 3S for k = 2(1)12 and 
n = 3(1)10, 12, 15, 20. Tables 5 and 6 give to 3S approximate upper 5% and 
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1% points for the ratio of maximum to minimum mean deviation for k = 
2(1)12 and m = 3(1)10, 12, 15, 30, 20, 60, ©. Tables 3 through 6 are based 
on the approximate distribution of Table 1 and the fact that the ratio of two 
independent x?’s has an F-distribution. Adequate indications of accuracy of 
the approximation are stated. 

W. J. Dixon 
University of Oregon 
Eugene, Oregon 


1238[K].—E. J. GumBeEL, Statistical Theory of Extreme Values and Some 
Practical Applications. NBS Applied Math. Series, no. 33, U. S. Gov. 
Printing Office, Washington, 1954, viii + 51 p., 20.1 X 25.9 cm. Price 
40 cents. 

This booklet contains four lectures given at the Bureau of Standards 
which are an excellent account of the theory of extreme values and their ap- 
plications by the leading expert in this field. For the exponential type of dis- 
tribution law for extreme values there is tabled (p. 29) for samples of N = 
15(5)50(10)100(50)300, 400, 500, 750, 1000, extremes, the expected mean of 
the reduced (standardized) extremes, the ratio of these to the population 
mean, their standard deviation, and the ratio of this to the population 
standard deviation, all to 5D. Interpolated values for the means and stand- 
ard deviations for N = 20(1)49 to 4D are also given (p. 31), these values 
being obtained from an empirically established linear relationship. 

aw & &. 


1239[K].—H. C. Hamaker, “ ‘Average confidence’ limits for binomial 
probabilities,’’ International Stat. Inst., Review, v. 21, 1953, p. 17-27. 
Given a random sample of m from a binomial population with true per- 

centage of defectives, p. If k items are defective, the usual upper and lower 

(1 — a) confidence limits for p are computed so that 


(1) Pr {pr’ <P < fe} [> 1, 


where the equality holds for only a few isolated cases and the confidence prob- 
ability is generally greater than 1 — a. The limits are computed from in- 
complete Beta or F-tables, so that the upper limit is 


(2) be = 1000, F/ (v2 + 01F), 


where v, = 2(k + 1) and ve = 2(m — k) are the degrees of freedom for F. 
A similar result holds for p’. 

This article shows how the confidence probability varies as p changes, 
for » = 25 and a = .20 (Figures 3A and 5B’). Another set of confidence 
limits is then presented for which the confidence probability is generally less 
than 1 — a (Figures 3B and 5A’). In this case the upper confidence limit is 
simply p,-1. 

The author proposes to use as average confidence limits: 


Pur = (Pe + Pe-1)/2; Pu = (pe! + Pe-1)/2 


(3) puo = 100 — pin = Pur/2 = (p1 — po)/4 
bi = 100 — pun = 0. 
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Tables of pu and p,, to 1D are given for m = 5(5)30, 40, 50, 75, 100(100) 
500; a = .02, .10, .20; k = 1(1)20, unless m < 20, when & goes to n. These 
confidence limits give confidence probabilities which do not fluctuate much 
from 1 — a; however, the confidence probability will sometimes be less than 
1—a. 

It should be mentioned that Smmon! presents a Bayesian method of 
averaging the probabilities; his limits are even narrower than those in the 
present article. It would be desirable to compare the true confidence prob- 
abilities for various p, using (2), (3) and the Simon charts. 

R. L. ANDERSON 
North Carolina State College 
Raleigh, North Carolina 


1L. E. Smuon, An Engineer's Manual of Statistical Methods. New York, 1941. 


1240[K].—N. L. Jonnson, “Some notes on the application of sequential 
tests in the analysis of variance,” Annals Math. Stat., v. 24, 1953, p. 
614-623. 


In the analysis of variance one often desires to test a linear hypothesis 
in some systematic model and such hypotheses are usually composite in 
nature. This paper considers the applicability of certain basic work by 
BARNARD! and Cox? in sequential tests for composite hypotheses to such 
cases. Such a sequential situation could be encountered if one treated the 
number of observations taken in each class as a random variable or if the 
number of classes studied in an experiment could be increased in a sequential 
fashion. In order to obtain a sequential test a restriction upon the manner of 
increasing the number of observations is needed. With such restrictions one 
can evolve a sequential test which is in a form similar to the usual sequential 
type, but the lack of adequate tables of required distributions makes it 
difficult to apply the theory in the general case. 

As a special case illustration the author considers a random model in 
one-way Classification, namely 


Xe = Ot ust 2 


in which the u’s are normal variables, each with 0 means and standard de- 
viations og. Here the hypothesis is concerned with the value of 6 = o*z/o’. 
If one is interested in a sequential test based on 6” > 8’ = 0, such a test 
exists, but to determine appropriate intervals for the tests, values of the 
criteria Gp and Gp are required. With the number 1 of groups of size k as the 
variable, Table Ia gives values of Ge and Ge to 3D for 6” = 1, 8’ = 0, 
a = B = .05 and k = 2(1)12, 15, 20, 30, 60, © and m = 2(1)12, 15, 20, 30, 
60. Table I, is similar but with a = 6B = .01. 
C. F. Kossack 

Purdue University 

Lafayette, Indiana 


1G. A. BARNARD, “The frequency justification of certain sequential tests,” Biometrika, 
v. 39, ee “yo 144-150. 
Ox 


v. 48, M980, p. 290-2 


uential tests for composite hypotheses,” Cambridge Philos. Soc., Proc. 
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1241[K].—Ranp Corp., Offset Circle Probabilities. Santa Monica, Calif., 
1953, 18 p., 16.5 X 24.1 cm. 


Consider a circle of radius rz with center located a distance D from the 
origin in the x,y-plane. Let g represent the probability that a sample value 
from the normal bivariate distribution with uy, = yw, = 0,0, = oy = o,p = 0 
falls outside this circle. 3D values of g are tabled as a function of D/e and 
(ra — D)/o. Here D/o = 0(.1)6(.5)10(1)20, © and (rz — D)/o =— 3.9 
(.1)4. This table was obtained from an unabridged 6D table of g which was 
computed jointly by the Institute for Nursrical Analysis, National Bureau 
of Standards, and the RAND Corporation. Card copies or printed listings of 
the unabridged table may be obtained by writing to the RAND Corporation. 


J. E. WatsH 
U. S. Naval Ordnance Test Station, Inyokern 


China Lake, California 


1242[K ].—S. Rosensau, ‘Tables for a nonparametric test of dispersion.” 
Annals Math. Stat., v. 24, 1953, p. 663-668. 


Given a sample of » values of a continuous random variable X, these 
tables give the upper 5% and 1% points for m,n = 1(1)50 of the distribution 
of the number of values in an additional sample of m values which lie outside 
the extreme values of the original sample. This distribution was obtained by 
WILks' in connection with the problem of two-sided distribution-free toler- 
ance limits. The suggested use in this paper is as a non-parametric test of 
dispersion, which places the emphasis on smaller values of m and n. How- 
ever, if neither the sample nor the unsampled portion of the population 
exceed 50, these tables can also be used to establish the level of tolerance 
limit provided by the extreme values of the sample. 

An approximation to the distribution for large and approximately equal 
values of m and 1 is also obtained. Since there is an apparent discrepancy 
between the 5% and 1% values obtained from this limiting distribution, and 
the values given by Wixks! for m = n = 100 (Table II, p. 406; 2 and N in 
WILKs notation), the reviewer computed exact values of the 5% and 1% 
points for m = n = 10(10)200(100)900. From these results it was noted (a) 
that the values given by WILKs are in error, and (b) that, while the 5% 
point of the limiting distribution is correct for m = m 2 25 (as can be seen 
from the tables), the exact 1% point does not reach that of the limiting 
distribution until at least m = n > 170. 

C. A. BENNETT 
General Electric Company 
Hanford Atomic Products Operation 
Richland, Washington 


1S. S. Wixks, “Statistical prediction with special reference to the problem of tolerance 
limits,” Ann. Math. Stat., v. 12, 1942, p. 400-409. 


1243[K ].—HERBERT SoLomon, “Distribution of the measure of a random 
two-dimensional set,’’ Annals Math. Stat., v. 24, 1953, p. 650-656. 
N(= 1,2) random circles of equal radius W are dropped according to a 

bivariate distribution with circular symmetry specified by the parameter o 

on a fixed circle of radius T with an aiming point at distance R from its 
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center. The measure of interest is the ratio of the covered area to the total 
area of the fixed circle. The author calculates for N = 1 the probability P¢ 
of getting at least C/W fraction coverage for specified values of T, R, N, a. 
Until now only the first two moments were known. 

Fig. 1 gives for N = 1 the contours of equal probability Pc = .05 
(.05).95 as a function of R for given values of W + aT, both in o units. 
The factor a is eliminated by the use of the second figure, which shows the 
relation between coverage 0.05 to 1.00 and the factor a(—.1 <a < .1) 
for fixed values of the quotient W/T = .5, .6, .75, .85, 1, 1.2, 2, 3, 5, 10, «. 
Thus the two figures represent the probability Pc of coverage C for given 
values of W, T and R (in o units). 

Four short tables give for N = 2 the lower and upper bounds Pc = 1 — 
(1 — Pc)? and Pe = P(c; + cz = C) to 3D. Table I holds for W/e = 1, 
T/o = 1, R=0 and C = .1(.1).9. Tables II and III deal with median 
coverage for W = T and W = 2T respectively. Table IV holds for W = 
.5T, C = .2. The graphs and tables demonstrate that a realistic decision 
can be made without resorting to random number devices. 

E. J. GUMBEL 
Columbia University 
New York, New York 


1244[K ].—G. Tacut1, “Tables of 5% and 1% points for the Polya-~-Eggen- 
berger’s distribution function,” Union of Japanese Scientists and Engi- 
neers, Reports of Statistical Application Research, v. 2, no. 1, 1952, p. 
27-32. 

Table I gives to 3S the minimum values of # such that 


F(k;h,d) = > (1/n!)h(h + d)---(h + (m — 1)d)(1 + d)-@rnole 


n=0 


is > a, where a = .95 or .99, h/d = .5(.5)15, 20, 30, 60, © and k = 1(1)25. 
Table II gives to 3S the maximum values of 4 such that 


F(k;h,d) 2 a 


where a = .95 or .99, d/h = 2., 1., .5(—.1).1, .05, 0 and k = 25(5)40(10)60, 
75, 100, 200, 500, ~. 
The method of preparing the tables and their accuracy is not given. 
I. R. SAVAGE 
National Bureau of Standards 
Statistical Engineering Laboratory 


1245[L ].—ApMIRALTY RESEARCH LaBoratTory, “Tables of f(g,a),”" A.R.L./ 
T.3/Maths. 2.7, 13 p., Teddington, Middlesex, England. 


4D tables of 


24 2 [T(m/2+1 
oo(- £5") [x00 +2 [Seren 


X u™?,F,(m/2; m + 1; — u) In(kg) cos ma} 
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for wu = 0,1, 0, k = .5(.5)3, a = 0(10°) 180°, g = 0(.2) go where go varies 
between 2 and 6 and is so chosen that the tabulation covers in each case the 
whole significant range of g. 


A. E. 


1246, L].—D. R. Bates, KATHLEEN LepsHAM & A. L. Stewart, ‘Wave 
functions of the hydrogen molecular ion,’’ Royal Soc. London, Phil. 
Trans., v. 246 A, 1953, p. 215-240. 

When the variables are separated in the Schrédinger equation for the 
electronic wave functions of H,* in elliptic coordinates, there result the two 
differential equations 








d dM m? 

(1) £ja- el + {-4 4 oa} no 
d dA m? 

(2) Kor oF] + {4+ 2m - pr —- bao. 


In the above, A is a separation constant, and 
(3) P? =— iRE, 


where R is the distance between the nuclei in atomic units, and E is the 
electronic energy in Rydbergs. 

Equation (1) is the familiar (oblate) spheroidal wave equation which has 
been studied.! Corresponding to a countable set of eigenvalues A (/,m,p) 
equation (1) has a solution of the form 


(4) M(m,p xu) = D fallsm,p) PM 4(u), 


where Pr+.(u) is the associated Legendre function and s is either always odd, 
or always even. Tables of the eigenvalues A and the coefficients f, have been 
published in! as functions of ». The authors give corresponding tables for 
the parameter R, which are presumably of direct interest to researchers in the 
field, and the tables are arranged in such a manner that the various orders of 
the eigenvalues are directly identified with the ‘‘united atom designations.” 
Ten such ‘‘states’”’ are given in the tabulation, namely 


nscg, nm =1,2,3; npou, m= 2,3,4; 3do,, Afou, 2pm., and 3dr. 


All ten correspond to parameters m = 0,1; / = 0,1,2,3. Corresponding to 
these states, the authors tabulate the functions p, ¢, —A’, and —E, for 
R = 0(.2)5(.5)9 or 10; 5D (Table 1), and the coefficients f, (Table 2). The 
latter are given to five significant figures for the dominant coefficients, and to 
six decimals (usually) for the others. According to the authors, these were 
obtained by interpolation from the tables in'. The parameters A’ and @ are 
defined below. 

The authors also tabulate eigenvalues and coefficients corresponding to 
equation (2). Following JAFFf,? the authors write 


(5) A(X) = (7 — 1)"?(\ + 1) exp (— pa)y(s) 








$$ rn, 
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with 
R 
ok Tee f= (A — 1)/A+1) 
A'’=A — P. 
The solution is written in the form 
(6) y= Des. 
t=0 


The solution exists corresponding to the eigenvalues of oc, for given p, m, 
and A’. In Table 3, the authors give the coefficients g, for the same range of 
R as in Tables 1 and 2; four decimals are given in the tabulation. The paper 
also contains contour diagrams of H,*, as well as some comparisons between 
the exact solution for certain of the wave functions and the “L.C.A.0” 
approximation. 

For theory and some tabulation relating to the general spheroidal wave 
function in the more accessible publications, see also references *-*:*. 

GERTRUDE BLANCH 

Consolidated Engineering Corp. 
Pasadena, California. 

1J. A. Srratron, P. M. Morse, L. J. Cau & R. A. Hutner, Elliptic Cylinder and 
Spheroidal Wave Functions. John Wiley and Sons, Inc., New York, 1941. 

2G. Jarré, Zeit. Phys., v. 87, 1934, p. 535 (author's citation). 

3A. Lerrner & R. D. Spence, ‘“‘The oblate spheroidal wave functions,” Franklin Inst., 
Jn., v. 249, no. 4, 1950, p. 299-321. 

4 J. MEIxNneR, Lamé’s Wave Functions of the Ellipsoid of Revolution, translated by Mary 
L. Manter, NACA, June, 1944. 


5C. J. Bouwxamp, “On spheroidal wave functions of order zero,” Jn. Math. Physics, 
v. 26, 1947, p. 79-92. 


1247[L].—S. CHANDRASEKHAR & DONNA ELBERT, “The roots of Y,(An) 
Jn(\) — Jn(An) Y,(A) = 0," Cambridge Phil. Soc., Proc., v. 50, 1954, 
p. 266-268. 
The authors tabulate to 5D the first root, A;, of the equation mentioned 
in the title, together with 7 to 9D values of 





Jn(A1), Y,(A1), Jn (dx), YnQum), oy Tn? (d1) >= 

for 7 = .2, m = 1(1)5; 9 = .3(.1).5, m = 1(1)6; 9 = .6, m = 1(1)8; and 
n = .8, nm = 1(1)12. For 7 = .5, m = 1(1)6 they also give the corresponding 
quantities for the second and third zeros. 


2 | Jn? (Arn) i} 


A. E. 


1248[L].—Cu. Vita Dunsk1, “Les fonctions de Bessel d’argument com- 
plexe xVj et les fonctions de Kelvin d’ordre zéro et 1,” Soc. Roy. Sci. 

Liége Bull., v. 23, 1954, p. 52-59. 

Table I gives ber x, bei x, ker x, kei x, ber; x, bei; x, ker; x, kei, x, and 
Table II gives the real and imaginary parts of Jo(xi#), Ho“ (xt), <4 J, (x4), 
4+ H,™ (xi#). Each of the two tables occupies one page. The entries have been 
computed, to various degrees of accuracy, from the asymptotic expansions 
for 25 values of x ranging from 10 to 72. 

A. E, 
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1249[L ].—Carit-Erik FrOpErG & Puitip Rasinowitz, “Tables of Cou- 
lomb wave functions,’”” NBS Report 3033, 1954, 3 + 40 p., 20.5 K 27cm. 


F and G being the regular and the irregular Coulomb functions as defined 
in RMT 982 [MTAC, v. 6, p. 92], in this report 


Ao(n) = [(1 — e**")/ (24m) }#, Ax(n) = 3(L* + 9°)-? Aral) 
Fy, ,(p) = Ax(n)p™ fr, a(0) 
Gr,4(e) = Ax(n)p** gz, 4(0). 
The report contains tables of f, f’, g, g’ and of the first five ‘‘reduced”’ deriva- 
1 oF 
tives Elan of these quantities for L = 0, 5, 7 = 1(1)10, p = 1(1)10. For 
other NBS tables of Coulomb wave functions see RMT 1091 [MTAC, v. 7, 
p. 101-102]. The present table was first announced as UMT 186 [MTAC, 
v. 8, p. 97], where further details are given. 
A. E. 


1250[L].—J. M. Hammerstey, “Tables of complete elliptic integrals,” 

NBS Jn. of Research, v. 50, 1953, p. 43. 

The functions tabulated are K, E, and 4x/K. The argument used is 
pb = 1/k. The tables are for p = 1(.01)1.3(.02)2; 10S. The last figure is 
doubtful. 

D. H. L. 


1251[L].—Atpo Mucc1a, “Sul calcolo dell’ integrale di Poisson,” Accad. 
Sci. Torino Cl. Sci. Fis. Mat. Nat., Atti, v. 87, 1953, p. 116-126. 
5D tables of 








4d t 
A(r) = 21n sin 5 : B(r) -f t cot = dt 
C(r) =— cot 5, D(r) =—- r cot 5 +2In sin 5 








for r =— mw, —3.12(.04) —.04(.01).04(.04)4.6. The tables were computed 
by Lia Errera of the computing center of the Politecnico di Torino. 
A. E. 


1252[L].—Uno Otsson, ‘‘Non-circular cylindrical gears,’’ Acta Polytechnica, 

Mechanical Engineering Series, v. 2, no. 10, Stockholm, 1953, p. 1-215. 

The author not only develops the theory for the construction of different 
types of non-circular gears but also furnishes directions for the production 
of the wheels. As a by-product of the investigation, he gives the following 
tables: 

Table I: The arc length of the hyperbola, 

Table II: The arc length of the ellipse with imaginary axes or the arc 
length of the sinh-curve, 

Table III: The arc length of the parabola, 

Table IV: The arc length of the exponential curve, 

Table V: The hyperbolic argument 4 as a function of the hyperbolic 
amplitude ¢. 











fe 


oeEEllC(<i rh CUCU TLC OCTUC~™ 





RECENT MATHEMATICAL TABLES 225 


Tables I and II are given for ¢ (the amplitude angle) = 0°(1°)90° and a 
(the eccentric angle) = 0°(5°)90°; Table III for g = 0°(1°)90° with A’, 
and Ao” + A;”’; Table IV for ¢ = 0°(.1°)10°, 0°(1°)90° and Table V for 
gy = 0°(.1°)90°. All tables are given for the most part for modified functions 
to about 5 significant figures. A short section describing the use of higher 
degree interpolation, necessary in the tables, is also given. 

IRENE A. STEGUN 
NBSCL 


1253[L].—L. J. Stater, “Some new results on equivalent products,” 
Cambridge Phil. Soc., Proc., v. 50, 1954, p. 394-403. 


Table 1 gives values, mostly to 8D or 8S, of 
| Ha - ae}: 
for a = — .9(.05).95 and g = 0(.05)1. 


Table 2 gives 8D or 8S values of the same quantity fora = 1,¢ = 0(.005) 
.89. 
The first table is stated to be accurate to 7D except where fewer figures 
are given, the second, to 8D. The tables were computed on the EDSAC. 
A. E. 


1254[L].—L. J. Stater, “The evaluation of the basic confluent hypergeo- 
metric functions,”’ Cambridge Phil. Soc., Proc., v. 50, 1954, p. 404-412. 


The functions discussed and tabulated in this paper are 


2. (9°). 9” 
(a) ZO. O- 
(2) «© (q*)n y" qin (at) 





n= (@)n (Q)n 


@). = (i -—@)(1 — g*)--- A — gt) 

Table 1 is of the function (2) for a = 0(.2)2, b = .2(.2)1,¢g= .9, y= 
(1 — g)x = .1(.1)1. 

Table 2 is of the same function, for the same values of a and 3, and for 
gq = .99, y = .01(.01).1. 

Table 3 is of the function (1) for the same values of a and 8}, and for 
q = .9, y = .1(.1).9. 

The values are mostly given to 7 or 8S: they were computed on the 
EDSAC. 


where 


A. E. 


1255[L ].—Rosert L. STERNBERG, “‘A general solution of the two-frequency 
modulation product problem. II. Tables of the functions Ama(h,k),” Jn. 
Math. Phys., v. 33, 1954, p. 68-79. 


“The purpose of this paper is to provide tables and evaluation methods 
for the functions 


Amn(h,k) = (2/x*) ff (cos u +k cos v — hk) cos mu du cos nv dv 
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for m,n = 0,1,|h| < 2,0 < k< 1 where R is the region defined by 
R:cosu+kcosu>h, O<uv<r 


and subsequently also to tabulate briefly the first six higher order functions 
Amn(h,k).” 

The functions arise in the problem indicated in the title and have been 
investigated in an earlier paper! where there are also further references. The 
tables are to 6D. 


A. E. 


1R. L. STERNBERG & H. KaurMan, “A general solution of the two-frequency modulation 
product problem. I,” Jn. Math. Phys., v. 32, 1953, p. 233-242. 


1256[L].—A. WaLTHER & H. UNGER, “‘Mathematische Zahlentafeln, nu- 
merische Untersuchung spezieller Funktionen,” Naturforschung und 
Medizin in Deutschland, 1939-1946, Band 3. Angewandte Mathematik, 
Teil I, p. 167-183. Verlag Chemie, Weinheim, 1953. 


This valuable report covers numerical tables of special functions com- 
puted in Germany during the period under review, and also analytical work 
of importance in connection with the numerical computation of special 
functions. The report consists of three parts: A. Introduction, B. Survey of 
developments and results, C. Bibliography. Each of parts B and C is divided 
in four sections: 1. elementary functions, 2. tables pertaining to astronomy 
and geodesy, 3. Bessel functions, 4. other higher transcendental functions. 
The bibliography lists some eighty items. Many of these have been published 
(and reviewed in MTAC). In the case of most of the tables which have not 
been published a brief description (in the style of the FMR Index) is 
appended. 


A. E. 


1257[L].—M. W. Wrrkes, “A table of Chapman’s grazing incidence 
integral Ch(x,x),’” Phys. Soc. Proc., v. 67 B, 1954, p. 304-309. 
This table gives 


x 
Ch(x,x) = x sin x f exp (x — x sin x/sin \) cosec? \ dX 
0 


to 3D for x = 50(50) 500(100)1000, x = 20°(1°) 100°, excluding values of 
Ch(x,x) which exceed 100. The table was prepared on EDSAC, using a 
Gauss 5-point formula repeatedly, with automatic adjustment of the length 
of the interval. A more complete discussion of the machine techniques used 
is to appear elsewhere. Various checks were applied: recalculation for 
x = 20°, 21°, 70°, differencing in both directions, and use of the relation 


Ch(x,x) + Ch(x,r — x) = 2 exp (x — xsin x) Ch(x sin x, $7). 


The values are expected to be correct to within a unit in the last place. 
5. F. 
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MATHEMATICAL TABLES—ERRATA 


240.—T. H. BLaKEsLey, A Table of Hyperbolic Cosines and Sines. London 
Physical Soc., 1890. 


“A comparison with a corrected version of a table by Hayasut has indi- 
cated 172 errors of one to three units of the last digit, as well as the 31 errors 
of 5 or more final units listed below. Errors in the argument are also included 
in the list. In addition 180 errors of between half and one final unit were 
found.” 


Page Function For Read 
1 arg 0.14 5 4 
arg 0.18 4 8 
cosh 0.43 .0938888 .0938833 
2 cosh 1.12 .7005670 .6955670 
sinh 1.12 .3642872 .3692872 
cosh 1.84 .2276799 .2276778 
3 cosh 2.03 .8727101 .8727109 
sinh 2.03 -7413746 -7413754 
cosh 2.04 .9103184 .9103190 
sinh 2.04 -7802896 -7802902 
cosh 2.05 .9483548 .9483180 
sinh 2.05 .8196198 .8195831 
cosh 2.06 .9867111 .9867119 
sinh 2.06 8592571 .8592579 
cosh 2.07 .0255038 .0255044 
sinh 2.07 .8993179 .8993187 
cosh 2.08 .0647395 .0646996 
sinh 2.08 .9398093 .9397694 
cosh 2.35 .2905196 .2904694 
sinh 2.35 .1951504 .1951003 
cosh 2.42 .6673910 .6673905 
sinh 2.42 .5784683 .5784688 
cosh 2.61 .8362940 .8362927 
sinh 2.61 -7627595 -7626582 
cosh 2.62 .9042644 .9042632 
sinh 2.62 .8314615 .8314604 
cosh 2.63 .9729254 .9729242 
sinh 2.63 -9008469 .9008457 
cosh 2.64 .0422838 .0422824 
sinh 2.64 .9709225 .9709212 
cosh 2.65 -1123463 .1123449 
sinh 2.65 .0416950 .0416937 
4 sinh 3.15 -6966062 .6466062 
arg 3.84 3 4 


Cara M. MunrorD 
The Univ. Math. Lab. 
Free School Lane 
Cambridge, England 
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241.—P. C. MAHALANosis, ‘Tables of random samples from a normal 

population,”’ Sankhga, v. 1, 1933, p. 303--328. 

We reconstructed MAHALANOBIS’ table of random normal deviates as 
follows: Using the random numbers from Tiprett’s table (Tracts for Com- 
puters, No. 15, Cambridge 1927) as probabilities, we read the corresponding, 
normal deviate from the Kelley Statistical Tables (1), New York, 1938. In 
doing so, we found 427 errors which are greater than 0.001, of which 219 
errors are greater than 0.01 with 132 errors greater than 0.1. In addition, the 
rows on ten pages are in the wrong order, and some blocks of numbers are 
totally incorrect, and much of page 306 is repeated on page 308. 

Plans are being made to publish revised tables in Sankhga. 


LORRAINE SCHWARTZ 
Statistical Laboratory 


University of California 
Berkeley, California. 


UNPUBLISHED MATHEMATICAL TABLES 


Unpublished tables of special functions are mentioned in RMT 1256 
and in Phil. Mag., s. 7, v. 45, 1954, p. 599-609. 


191[A].—F. L. Mrxsa, A Table of Binomial Coefficients for N =1 to 
N = 100. Typewritten manuscript of 41 p. deposited in the UMT 
FILE. Other copies are obtainable gratis from the author. 

This table gives exact values of the binomial coefficients for the first 
hundred integer powers. It is thus equivalent to RMT 1234. A comparison 
of the two tables has not been made. 

F. L. MIKsa 

613 Spring St. 

Aurora, Illinois. 


192[F].—-A. Ferrier, Etude de 5n? + 1 et 5n‘ + 1. 13 typewritten pages 
deposited in the UMT FILE. 
This study is in three parts. The first part gives the solutions of the 
congruences 
5n*? + 1 =0 (mod 9), p < 12000 
Sn* + 1=0 (mod p), p < 12000 


The second part is a list of numbers 5n? + 1 for m < 5000 which havea prime 
factor exceeding 10’. The last part is a factor table of 5”* + 1 complete to 
n = 100 and partially complete tom = 1000. 

A. FERRIER 
College de Cusset 
Allier, France 


193[F].—F. L. Mixsa, Table of All Primitive Pythagorean Triangles Having 
Equal Areas Below 97017 29310. Twelve mimeographed leaves deposited 
in the UMT FILE. 


Fermat noted that if m and m are integers then the two triangles whose 
sides are equal to the absolute values of 
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(3n — m)*n? + m?(m + n)?, (3n — m)*n? — m?(m + n)?, 
2mn(2mn — m? + 3n?) 
and 


(3n + m)*n? + m?(m — n)*?, (3n + m)*n? — m?(m — n)?, 
2mn(2mn + m? — 3n?*) 


constitute a pair of Pythagorean triangles of equal area. Not all such pairs 
are given by this rule, however. The table gives all pairs of such triangles 
which are primitive whose areas do not exceed 10”, arranged according to 
increasing area. Besides the area, which is given also in factored form, the 
three sides and the Pythagorean generators of the triangles are given. There 
are 149 such pairs of triangles and one triplet. Of these only 81 are given by 
Fermat’s rule. The triplet is made up of triangles whose sides are 


(8580,3059,9109) (1380,19019,19069) (5852,4485,7373) 


Each has an area of 13123110. 


F. L. MrKsa 
613 Spring St. 
Aurora, Illinois 


194[F].—F. L. Mixsa, Counts of Primitive Pythagorean Triangles by Area. 
7 Hectographed leaves deposited in the UMT FILE. 


This list gives the number of primitive Pythagorean triangles whose areas 
lie between various limits, up to 10". More precisely the intervals are 


n-10' <A < (n + 1)-108 


for 
k=6, n=1(1)9 
k=7, n= 1(1)99 
k= 8, n = 10(1)99 


The total count for triangles with area less than 10" is 52482. 


F. L. MrKsa 
613 Spring St. 
Aurora, Illinois 


195[G].—R. L. Bivins, N. MEtropo.is, P. R. Stein & M. B. WELLs, 
Character Tables for Symmetric Groups of Degree 15 and 16. One reel of 
microfilm deposited in the UMT FILE. 

This table is the one referred to in the authors’ paper MTAC, v. 8, p. 

212-216. 


196[K ].—Mary G. NatreLta, Critical Values for Wilcoxon Rank Sum 
Test. Six mimeographed leaves deposited in the UMT FILe. 
Critical values of smaller rank sum for the Wilcoxon test for the 5,2, 
and 1 percent significance levels are given. These are compiled from various 
published tables with corrections of all known errata in the originals. 


C. EISENHART 
NBSSEC 
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197[L].—A. OpLer, Table Related io the Error Function. 20 leaves tabulated 
from punched cards deposited in the UMT FILE. 


The table gives 6S values of the integral 


2 exp (2?) f exp (—?*) dt 


for z = — 4.5(.1)100, together with first differences. The table is probably 
good to only 5S for most of the range. 

A. OPLER 
Dow Chemica! Co. 
Pittsburg, California 


AUTOMATIC COMPUTING MACHINERY 


Edited by the Staff of the National Applied Mathematics Laboratory of the National 
Bureau of Standards. Correspondence regarding the Section should be directed to Dr. E. W. 
Cannon, 415 South Building, National Bureau of Standards, Washington 25, D. C. 


BIBLIOGRAPHY OF CODING PROCEDURE 


68. J. L. Jones, A Survey of Automatic Coding Techniques for Digital Com- 
puters. 

This MIT master’s thesis is one of the first publications on the gen- 
eral subject of automatic coding. It includes chapters on the basic forms 
of automatic coding, a survey of present techniques, and a survey of 
contemplated techniques with comments. It includes an appendix of 
additional applications of the interpretive method and an appendix 
listing computer groups interested in automatic coding. 

69. C. C. Etcot, U. S. Naval Ordnance Laboratory, White Oak, Md., 
Single v. Triple Address Computing Machines (NAVORD Report 2741). 

By utilizing a slightly idealized notion of a single address computer, 
the author obtains a partial answer to the question of whether a single 
address or a triple address computing machine requires fewer words to 
specify a sequence of instructions. While he finds that fewer words are 
required to code by means of an “‘idealized’’ single address machine, it 
is interesting to note that the pseudo codes for single address machines 
tend to be three-address codes. 

70. H. B. Curry, U. S. Naval Ordnance Laboratory, White Oak, Md. A 
Program Composition Technique as Applied to Inverse Interpolation 
(NOL Memorandum 10337). 

71. J. W. Backus, International Business Machines Corporation. Opera- 
tional Characteristics of the IBM 704. 

The IBM 704 differs from the 701 in that among other things four 
floating point operations have been added along with automatic address 
modification. Here floating binary point representation similar to that 
presently used in pseudo-codes can be utilized directly by the computer. 

72. W. Brown, J. DETurK, H. Harner, E. Lewis, University of Michigan. 
The MIDSAC Computer. 
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The MIDSAC is a high-speed digital computer designed for use in an 
automatic control system. This 55 page booklet includes chapters on 
programming rules and operational characteristics, a detailed description 
of the computer, and the physical construction of the computer. 

73. G. E. ForsyTtHE, SWAC computes all 126 distinct semigroups of order 4. 
Manuscript 5 p. 

The problem of finding all (non-isomorphic) groups of a given order 
is known to be very difficult unless the order is of a very special nature. 
The possibility of using electronic computing machines in this connection 
has occurred to a number of people, but until now no results have ap- 
peared. An idea of the complexity of the problem is obtained from a 
study of this computation of all distinct semigroups of order 4. By semi- 
group is meant a set of elements closed under an associative composition 
law, called multiplication for simplicity; two such sets are distinct if 
they are neither isomorphic nor anti-isomorphic. (Two systems are 
called isomorphic (anti-isomorphic) if a one-to-one mapping X — X’ 
exists such that (AB)’ = A’B’ (= B’A’).) 

For machine purposes the elements of the semigroups are denoted 
by 0,1,2,3. The semigroup is then completely described by its multi- 
plication table, i.e., the matrix (a;;), 7,7 = 0,1,2,3, where a; is an integer 
which specifies the product of the elements 7,7. Four SWAC codes have 
been prepared by Forsythe in order to make a table on punched cards 
of all 3492 semigroups of order 4 and to pick out the 126 distinct among 
them (in previous work only 121 have been found). The 4 rows of the 
multiplication table are listed in one line consecutively, e.g. 0000 0000 
0023 0032. The distinct semigroups are represented by normal forms dis- 
tinguished through lexicographic order. 

The first code checks on the associativity. A direct test of all the 4'* 
multiplication tables would take SWAC about 13 years. However, by 
testing partial tables, and then rejecting blocks, the time was cut down 
to less than 80 minutes. The second code is for obtaining the normal form. 
The third code picks out the distinct semigroups. The fourth code con- 

verts the 126 tables from base 2 notation to base 4. 

The codes used are on file at UCNAR and NBSCL. 

Oxtca TAUSSKY 
NBS 


BIBLIOGRAPHY Z 


1147. H. J. Gray, Jr., “Numerical methods in real-time simulation,” 
Quart. of Applied Math., v. 12, 1954, p. 133-140. 


Under the sponsorship of the Office of Naval Research, Special Devices 
Center, the University of Pennsylvania conducted a study to determine the 
feasibility of building a digital real-time simulator for a system described by 
thirteen first-order non-linear differential equations. It was found necessary 
to improve either the speed of the digital machine or the quality of the 
method of solution. The paper deals with the latter problem and especially 
with the buildup of errors. 

The paper is introductory in nature, and only linear constant-parameter 
systems of ordinary differential equations are considered. The characteristic 
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equation associated with the corresponding difference equation is found and 
solved by conformal mapping. A diagram is given showing the method of 
solution in a case of practical interest and not requiring high accuracy. 

E. W. C. 


1148. Aston S. HousEHOLDER, “‘Generation of errors in digital computa- 

tions,” Amer. Math. Soc., Bulletin, v. 60, 1954, p. 234-247. 

This paper is an address delivered before the Spartanburg meeting of the 
Society on November 28, 1953. 

The problem of computing the numerical values of f(x) is considered. 
In general one must use, in the computation, an approximation x* to x, and 
an approximation f,(x*) to the function f(x*). The error is expressed as a 
sum of three independent components: f(x) — f(x*), f(x*) — fa(x*) and 
fa(x*) — f*(x*), called, respectively, the propagated error, the residual error 
and the generated error. The article treats the generated error. It is stated 
that differences among high-speed computing machinery are such that each 
machine demands a theory on the generated error. A specific computing 
system, the ORACLE, is described and the generated error for certain com- 
putational sequences is discussed. 


E. W. C. 


1149. L. Knicut, “Valve reliability in digital calculating machines,” 

Electronic Engineering, January, 1954, p. 9-13. 

A survey is presented of some considerations in selecting vacuum tube 
types, designing circuits to provide useful tolerences for the tubes, testing 
them, and maintaining vacuum tube equipment in order to achieve a level 
of reliable operation demanded by digital calculating machines. Special con- 
sideration is given to effects that develop as tubes age. Results of some 
vacuum tube experience with the Hollerith type 541 multiplier and the 
Manchester University computer are given. 


RussELL A. KirscH 
NBSCL 


1150. T. PEARCEy, Use of Punched Cards for Fourier Synthesis. Division of 
Radiophysics, Commonwealth Scientific and Industrial Research Organ- 
ization, RPR 121. Melbourne, November 1953, p. 1-20. 

The synthesizing of series of orthogonal functions is described in detail 
for the case of two-dimensional Fourier synthesis. Hand selection and refiling 
of cards from two master files of about 4000 cards are performed. The master 
files should allow for adequate spacing of the variables and positive as well 
as negative coefficients. The accuracy obtainable by the described method is 
three digits. The only equipment necessary is a tabulator for summing, a 
sorter and possibly a collator. 


STANLEY PRUSCH 
NBSCL 


1151. R. M. Saunpers, “Digital computers as an aid in electric machine 
design,”’ Communications & Electronics, May 1954, p. 189. 
An S.O.S. from Industry to the Electronic Brain to take over the drudg- 
ery of re-iterative computations required for designing electric machines. 
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There are tables showing the variables in the design of synchronous machines 
and a table of operations and sample calculations for ascertaining perform- 
ance of a 3-hp. induction motor. The “digital approach” and its advantage 
over hand computations is discussed. 

Discussion, pro and con, by Ropert L. FrttmoreE and C. G. VEINoTT is 
followed with a rebuttal by Saunders. 

A. R. Cock 
NBSCL 
NEws 


Numerical Analysis Research, UCLA.—Numerical Analysis Research, UCLA, is con- 
tinuing much of the research program of the Institute for Numerical Analysis (INA) of the 
National Bureau of Standards; INA was disbanded on 30 June 1954 [MTAC, v. 8, p. 178]. 

The new organization, administered by the UCLA Mathematics Department, will re- 
ceive support from the Office of Naval Research and the Office of Ordnance Research. Its 
primary mission includes pertinent fundamental research in mathematics and science and 
research in the application of computers to problems occurring in science and other applied 
fields. It will be aided in attaining these objectives by the use of the National Bureau of 
Standards Western Automatic Computer (SWAC) and other equipment which has been 
loaned to UCLA. In addition, the new organization has the use of INA’s library. Most of the 
research staff of INA has joined the new organization. 

The organization offers a training program in the efficient application of high-speed 
digital computers. A number of graduate assistantships offered in cooperation with various 
university departments are available each year. Seminars, and formal and informal courses 
in numerical analysis are conducted. 


Correspondence should be addressed to: 


Numerical Analysis Research 
University of California 
405 Hilgard Avenue 
Los Angeles 24, California 


OTHER AIDS TO COMPUTATION 


A Useful Technique in Programming for 
Analog Computers 


In programming for analog computers programmers have found certain 
types of linear differential equations difficult to handle even with the use of 
arbitrary function generators. In this paper a technique is presented which 
enables the analog computer to solve a heretofore difficult set of problems. 
The technique, while most useful in linear problems, is also useful in other 
problems. In many cases the technique will reduce the sensitivity of the 
solution to noise. 

Let us consider the general ordinary differential equation 


(1) Sony’, . + x) =0 
with the parametric transformation 


(2) $(x,z,2, ia t) = 0. 
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One sets up feedback circuits to generate x and functions of x and uses 
these to generate y. In the case ¢ = x — ¢ = 0 one gets the usual setup. 
The special case ¢ = x + e** = 0 has several advantages in certain cases. 
Two examples following show its use: 

Consider the equation 


(3) (2x8 + x)y’ + (3x71 + x-”)y = 2x-48, 
For x > 1, let x = e*. Then one gets 
(4) y = [ 2e-45¢ ote y (3e7-1# a e™)/a a 2e-#/8) 


which is easy to set up, generating the exponentials in individual integrator- 
pot-amplifier loops. Of course, one plots y against e‘, not ¢. 

An example for which this technique is useful in a non-linear equation is 
the following : 


(5) Ay” + By’ + Cyle?' = 0. 
Let v = dy/dx and let y = e*. Then dx/dt = e*/v and thus (5) becomes 
(6) » = ae* + Bu exp (3¢ — De-‘), 
where 
a=—B/A, B=-—C/A, 


which can be generated without the use of special function generators. 

A parametric substitution as in the above examples can frequently 
eliminate or considerably reduce the need for special function generators. 
Other parameterizations can be used, but that which is best can only be 
determined by an investigation of the specific problem under consideration. 


J. E. MAXFIELD 
Naval Ordnance Test Station 
China Lake, California 


BIBLIOGRAPHY Z 


1152. S. Amitsur, U. OPPENHEIM & A. Many, “An electrical computer for 
the solution of linear simultaneous equations,’’ Rev. Sci. Instruments, 

v. 24, 1953, p. 112-116. 

A variation of the Many-Meiboom network! is given which will solve a 
symmetrical system of linear equations. The given system of equations is 
analogous to the statement of Kirchhoff’s laws for a circuit of pure react- 
ances. Transformers permit the representation of negative nondiagonal 
quantities by condensers but diagonal negative quantities necessitate the 
use of coils. The problem of phase difference between the input and output 
alternating voltages is discussed in detail. 


F. J. M. 


1A. Many & S. MErsoow, “An electrical network for determining the eigenvalues and 
eigenvectors of a real symmetric matrix,”’ Rev. Sci. Instruments, v. 18, 1947, p. 831-836. 
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1153. V. C. ANDERSON & Partie Rupnick, “A thermistor bridge correlator,” 
Rev. Sci. Instruments, v. 24, 1953, p. 360-361. 
Thermistors are used to obtain the squares of voltages for use in a 
quarter square multiplier. 


F. J. M. 


1154. N. F. Barper, “Harmonic analysis,” Rev. Sct. Instruments, v. 24, 
1953, p. 329. 


An improvement in the method of F. B. DANTELs! is suggested. 


1F. B. Danrets, “A planimetric method of harmonic analysis,” Rev. Sci. Instruments, 
v. 23, 1952, p. 369-370; MTAC, v. 7, p. 120. 


1155. M. L. Boas, “‘Nomograms used in analysis of data in the reaction 
Li®(”,a)H*,”’ Rev. Sct. Instruments, v. 24, 1953, p. 356-359. 


1156. S. L. Brown & H. D. Scowetman, “The application of the Laplace 
transformation and a mechanical harmonic synthesizer in the analysis of 
electric circuits,” Rev. Sci. Instruments, v. 24, 1953, p. 375-379. 

The harmonic synthesizer is a purely mechanical one, using gear ratios, 
scotch yokes and a chain adder. The harmonic synthesizer is used to plot 
the values of P(z) on a circle |z| = r. Thus the real part of P(z) is given by 


> a:r* cos ky and the imaginary part by a similar expression. The synthe- 
k=0 

sizer is used to plot the complex number P(z) and the number of times this 
graph loops the origin is the number of roots of P(z) within the circle. Two 
examples are given with detailed references to specific circuits. 


F. J. M. 


1157. G. L. BROWNELL, R. V. Caviccut & K. E. Perry, “An electrical 
analog for analysis of compartmental biological systems,” Rev. Sct. 
Instruments, v. 24, 1953, p. 704-710. 

An analogy is set up between a biological system consisting of fluids in 
various compartments and an electrical circuit consisting of condensers 
connected by resistors. The condensers receive an initial charge at intervals 
of a tenth of a second and the transient response of the system is indicated 
on an oscilloscope. The system is used to analyze the results of tracer experi- 
ments. Models or simply diffusion constant values are hypothecated and 
tested by the analyzer and the results compared with experimental values 
obtained by the use of traces of radioactive elements. 


F. J. M. 


1158. A. B. CampeL, & Paut J. SCHNEIDER, ‘Membrane apparatus for 
analogic experiments,”’ Rev. Sct. Instruments, v. 24, 1953, p. 513-514. 
A soap film apparatus for solving the Laplacian equation is described. 


1159. I. CepERBAuM, ‘“‘Note on high speed product integrator,”’ Rev. Sci. 
Instruments, v. 24, 1953, p. 1072-1073. 
The author obtains an exact expression for the error due to the variation 
of y in the Macnee set up described in MACNEE (1163) in the case in which 
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the kernel K is trigonometric and y varies either continuously or in small 
discrete steps. 


F. J. M. 


1160. C. B. Crump, Jr., “Engineering uses of analog computing machines,”’ 

Mechanical Engineering, v. 74, 1952, p. 635-639. 

This is an introductory article which points out the nature of differential 
analyzers, both mechanical and electrical, various possible applications and 
a brief description of how such a computer is used. 

F. J. M. 


1161. C. A. LuDEKE & R. T. Evans, “A coupling analog for nonlinear svs- 
tems with more than one degree of freedom,” Jn. Appl. Physics, v. 24, 
1953, p. 119-122. 

C. A. LupEKE & C. L. Morrison, “Analog computer elements for solving 
nonlinear differential equations,” Jn. Appl. Physics, v. 24, 1953, p. 
243-248. 

The objective of the work of these papers is to develop a differential 
analyzer, whose elementary component contains a representation of a non- 
linear differential equation of the second order, for instance 


a0 
an +) = 0. 


In the second paper, which seems to precede the first logically, such an ele- 
ment is developed in the form of a coil suspended as a pendulum within an- 
other coil. The angle of rotation of the coil is the dependent variable and it 
can be used in turn to rotate a mask controlling the amount of light which 
falls on a photo tube. The output of the latter is amplified and the resulting 
current passes through the coil on the pendulum. This yields a possible 
representation of f(@) as a restoring force dependent upon the dependent 
variable, 6. A number of examples of periodic motions are quoted and good 
accuracy in the frequency is indicated. In the first paper, coupling between 
two such units by mounting two coils on each pendulum is described. 
F. J. M. 


1162. F. J. McDonat, “A Fourier analyzer,’ Rev. Sct. Instruments, v. 24, 

1953, p. 272-276. 

This instrument was developed to provide frequency analysis of transient 
signals, obtained on a multichannel recorder with band width of 5 to 200 cps, 
during seismic prospecting for oil. The record on variable area film is read by 
a photoelectric cell, and thus the function to be transformed is obtained as a 
voltage. A scotch yoke produces the sines and cosines in a form which yields 
a shaft rotation and hence a linear potentiometer can be used as a multiplier. 
Integration is by a tachometer type servo integrator and the output appears 
on a revolutions counter. 


F. J. M. 


1163. A. B. Macneg, “A high speed product integrator,”’ Rev. Sci. Instru- 
ments, v. 24, 1953, p. 207-211, 
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The integrals 


T v 
F(y) = f fOKOMd or f fOKOAdt 


are evaluated by means of high speed analog computer elements for various 
values of y. The variable y is given a sequence of values and by a suitable 
switching procedure the values of F(y) are evaluated at a rate of 60 per 
second and plotted with a reference value on a cathode ray tube. The func- 
tions f(#) and K(y,t) for y fixed are obtained by function generators or by 
differential analyzer techniques. During the generation of F(y), y which 
should remain fixed does actually vary. The error due to this cause and the 
errors due to the multiplier and functior generator are discussed. 
F. J. M. 


1164. J. R. Racazzin1t & G. REYNOLDs, “The electronic complex plane 

scanner,” Rev. Sci. Instruments, v. 24, 1953, p. 523-527. 

A rational function F(z) is supposed given in the form of a quotient of 
products of linear factors z — ),. It is required to obtain log | F(jw)| and 
arg F(jw) where j = V— 1 and w is real. Log | F(jw)| is the sum of + log 
| (jw — d,)|. The quantity jw — i; is represented as a complex voltage and 
log | (jw — ,)| is obtained by a suitable circuit. The argument of F(jw) is 
obtained by differentiation using the Cauchy-Riemann equations. 

F. J. M. 


1165. J. TADAYON, ‘“‘Measurement of the angle between two curves,” Reo. 
Sci. Instruments, v. 24, 1953, p. 871-872. 
An optical instrument for measuring the angle between two curves on a 
graph is based on the mirror principle for finding the normal to a curve. 


F. J. M. 


NOTES 


168.—A PRACTICAL REFUTATION OF THE ITERATION METHOD FOR THE 
ALGEBRAIC EIGENPROBLEM. In the second part of my paper on algebraic 
eigenproblems? I have proved that the computation by means of the forma- 
tion of the characteristic equation requires less computational work than the 
iteration method, and that this holds even when nothing but the first eigen- 
value has to be calculated. This advantage grows with every accessory eigen- 
value or vector. Further one has no trouble with deflation which requires a 
lot of multiplications. Also one can compute every eigenvalue and vector 
apart from the others, and do this to any desired accuracy, by the more 
powerful methods for algebraic equations. At last, there arise no difficulties 
from an unfavorable quotient of the two dominant eigenvalues. 

The reason these facts are not yet universally acknowledged is that the 
iteration method seems to be simpler and more mechanical in its application. 
But one has to consider that iteration does not converge quickly enough in 
practice, unless the quotient of the two dominant eigenvalues is } or less. 
This last will be rarely the case. For the eigenvalues must lie somewhere 
between two circles around the origin in the complex plane. The radius of 
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the inner circle has a value quite different from zero (if the determinant of 
the matrix is ~ 0). In this ring can lie only some of the m eigenvalues having 
a quotient less than 4. The other eigenvalues must cluster around some 
circles. In proportion as the situation is more favourable for e.g. the two 
dominant eigenvalues, the smaller the room will be for the other ones. There 
is no escape from this situation, and sooner or later the then dominant 
eigenvalues will lie close together. 

That iteration may give trouble in the most simple cases is a well known 
fact from numerous calculations. Yet the following example of a symmetric 
4 X 4 matrix, the eigenproblem of which is therefore real, is amazing. 


2 1 3 4 
1 =—3 1 5 
3 1 6 =—2 
4 5s -2 —1 


A= 


Let » be the dominant eigenvalue. Following AITKEN,! we begin the iteration 
with the usual starting vector v = (1,1,1,1)’ and so compute vm, = Avm_i. 
The quotient uw» of the first components of v, and vm_1 yields for u the ‘‘ap- 
proximations’”’ : 
a Be = (ha ae = ‘ae Me 7.85, Msp = 7.56, Me = 8.0, ha = 7.6, 
us = 8.08, wo = 7.667, wio = 8.122, wir = 7.683. 

The convergence is not impressive, especially if one considers that the se- 
quence S converges to about —8 (minus eight)! On the contrary, S seems to 
diverge or to “converge” to two limits. 

To analyze the sequence S more deeply, all eigenvalues and vectors were 
computed from the characteristic equation, the vectors by the method de- 
scribed in BopEwIG,? part II, section 3, p. 1-3. These are to 8D 


wp = —8.02857835, nu =7.93290472, uv = 5.66886436, wu = —1.57319074. 
x® = (1, 2.50146030, —0.75773064, —2.56421169)’ 

x = (1, 0.37781815, 1.38662122, 0.34880573)’ 

x*® = (1, 0.95700152, —1.42046826, 1.74331693)’ 

x® = (1, —0.90709211, —0.37759122, —0.38331238)’. 


Since S does converge to —8.0...the situation seems to be that the yu; 
diverge in the beginning. The 2; decrease and the yuo;+; increase, and in order 
to converge to —8, the ue; go through zero and the uo; through ©. By a 
special method the index has been determined where the signs change. This 
is the case at the 363rd or 364th iteration. 

To demonstrate this conclusion we have computed the values of yu; which 
would have resulted if the iterations had been really effected. The 100th to 
102nd iteration would yield the ‘‘approximations’”’ : 


Hi00 = 7.28514, M101 = 8.64677. 
The 200th to 202nd iteration would yield: 
Hoo0 = 5.96936, peor = 10.57380. 
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The 300th to 302nd would be: 
200 = 2.86532, sor = 22.13220. 


Our above conjecture concerning the y2; and 2,41 is therefore confirmed, as is 
also apparent from theoretical reasons. Near the critical point 363 we should 
get 


bso = 0.13587, uss1 = 468.66407, use2 = 0.04022, uoss = 1583.336, 
Moe = — 0.055448, wos = — 1148.735. 


The 400th to 402nd iteration would give 
Maoo = — 1.75101, sor = — 36.46902, 


the 800th to 802nd: 
vsoo = — 7.94348. ps0 — 8.11356, 


and at last the 1200th to 1202nd: 
#1200 = — 8.02787279, i201 = — 8.02927758. 


Thus 1200 iterations will scarcely yield 4 figures of the dominant eigenvalue! 
And this for a simple matrix of order 4. 

This amazing conduct can afterwards be explained from the knowledge 
of the eigenvalues and vectors. The slow convergence has two sources. Not 
only have the dominant eigenvalues nearly the same absolute value, but 
also the starting vector v is nearly orthogonal to x™. In fact writing x™ in 
the approximate form (4,10,—3,—10)’, the cosine between v and x" is only 
1/30 which is about cos 88°. 

Yet the situation would not considerably improve if another starting 
vector, e.g. (1,0,0,0)’ would be taken. Nor should we waste time by com- 
puting the eigenproblem of matrices of the form A + ¢ with varying c’s. 
The time is better spent by computing the characteristic equation as is 
pointed out in BopEwic? (Part II in Zusammenfassung”). Even the use of 
quadratic equations would not be very efficient in our case as »™ lies close 
to n». Nor would this be the case for other matrices if one wants to have 
more than the two dominant eigenvalues, and even then the determination 
of the vectors is far from agreeable. 

We owe to von Miszs? the discovery of the iteration method for finding 
the dominant value and vector and to AITKEN its application to all compli- 
cated cases more or less, and to HOTELLING the discovery of deflation. But 
times have changed. In our days also the higher eigenvalues and vectors 
must be computed and with extra accuracy. A deeper analysis shows that 
the iteration method in its present form is not appropriate and that even 
in cases when, by exception, the first eigenvalue is furnished quickly the later 
eigenvalues present the greater troubles. 

So the whole eigenproblem must be considered anew. 

E. BoDEWIG 
Borweg 7 
The Hague 
Holland 
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1A. C. AITKEN, “The evaluation of the latent roots and latent vectors of a matrix,” 
Roy. Soc. Edinburgh, Proc., v. 57, 1936-1937, p. 269-304. 
*E. Bopewic, “Bericht iiber die Methoden zur numerischen Lésung algebraischer 
Eigenwertprobleme,” Seminario matematico e fisico dell’ Universita di Modena, Atti, v. 4, 
1949-1950, v. 5, 1950-1951. 
+R. von Miszs, “Praktische Methoden zur Gleichungsauflésung,” Zeit. angw. Math. 
Mech., v. 9, 1929, p. 62-77. 
169.—ANALYTICAL APPROXIMATIONS. [See also NoTE 157. ] The follow- 
ing are approximations for the exponential integral and certain Bessel 
functions. a 
(56) O0<x«<i1 |e] max = 2 X 107 
— Ei(—x) + log. x = — .57721566 + .99999193 x 
— .24991055 x? + .05519968 x* — .00976004 x* + .00107857 x5 or! 
(57) —3<x<3 |elmax = 1077 a 
Jo(x) = 1 — 2.2499997 (x/3)? + 1.2656208 (x/3)* pore 
— .3163866 (x/3)® + .0444479 (x/3)® — .0039444 (x/3)™ on 
+ .0002100 (x/3)” Si; 
(58) O< x<3 |e|max = 2 X 10-* lik 
¥,(x) —2 logs = Jo(+e) + .36746691 + .60559366 (x/3)* af 
T 
— .74350384 (x/3)* + .25300117 (x/3)* — .04261214 (x/3)® sil 
+ .00427916 (x/3)* — .00024846 (x/3)" te: 
(59,60) 3<x< ™ 
Jo(x) = x~* fo(3/x) cos {x — go(3/x)} ex 
Yo(x) = x~* fo(3/x) sin {x — go(3/x)} va 
|€| max = 10-* ro 
fo(3/x) = .79788456 — .00000077 (3/x) — .00552740 (3/x)? pr 
— .00009512 (3/x)* + .00137237 (3/x)* — .00072805 (3/x)§ w 
+ .00014476 (3/x)*® N 
|e] max = 5 X 10-8 : sp 
¢o(3/x) = .78539816 + .04166397 (3/x) + .00003954 (3/x)? st 
— .00262573 (3/x)* + .00054125 (3/x)* + .00029333 (3/x)® N 
— .00013558 (3/x)® 
(61) 0<x<3 |e] max = 5 X 10° 








Ji(x)/x = 5 — .56249985 (x/3)? + .21093573 (x/3)4 
— .03954289 (x/3)® + .00443319 (x/3)? — .00031761 (x/3)™ 
+ .00001109 (x/3)” 


(62) O<x<3 |elmx = 5 X 10-8 


Yi(x) — = hog. re) = — 6366198 + .2212091 (x/3)? 
+ 2.1682709 (x/3)* — 1.3164827 (x/3)* + .3123951 (x/3) 
— .0400976 (x/3)© + .0027873 (x/3)" 

(63, 64) 3<x< © 
Ji(x) = x-* f1(3/x) sin {x — ¢(3/x)} 


Yi(x) = — x4 f,(3/x) cos {x — ¢1(3/x)} 
|e] max = 3 X 10-8 


A 
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f1(3/x) = .79788456 + .00000156 (3/x) + .01659667 (3/x)? 
+ .00017105 (3/x)* — .00249511 (3/x)* + .00113653 (3/x)® 
— .00020033 (3/x)® 
|e] mex = 8 X 10-* 
¢1(3/x) = .78539816 — .12499612 (3/x) — .00005650 (3/x)? 
+ .00637879 (3/x)* — .00074348 (3/x)* — .00079824 (3/x)® 
+ .00029166 (3/x)® 
E. E. ALLEN 
Shell Development Co. 
Houston, 25, Texas 


170.—A SIEVE PROBLEM ON “PsEUDO-SQUARES.” The following problem 
originated by KRaITCHIK,' and extended by LEHMER? by a special sieve, has 
recently been further extended by the SWAC. Let p be a prime. Let N, be 
the least positive non-square integer of the form 8x + 1 that is a quadratic 
residue of all primes < ». In this definition, zero is not counted as a quad- 
ratic residue so that N, is not allowed to be divisible by any primes < p. 
Since squares are quadratic residues of any prime, the numbers N, behave 
like squares and may be called pseudo-squares. This fact makes the problem 
of discovering pseudo-squares not only a sifting problem but also one of re- 
jecting squares. Thus the problem is unsuitable for a high speed sieve alone 
since the output would be mostly squares, each of which would have to be 
tested by a more elaborate arithmetic unit. The problem is therefore one for 
an all-purpose computer programmed for sifting.* 

Since for p > 3, N, must be of the form 24x + 1. One may proceed to 
exclude values of x using prime moduli between 5 and # inclusive. For every 
value of x not excluded the machine is programmed to extract the square 
root of 24x + 1. If this is a perfect square, the machine returns to the sifting 
program for the next value of x. Fortunately the early part of the program, 
where the squares come thick and fast, had already been carried? as far as 
Nei = 48473881 in 1928 so that when programmed for the SWAC the routine 
spends most of its time sifting. Actually, for the record, the SWAC was in- 
structed to print out every 64th square it produced. The complete table of 
N, for p < 83 is as follows. 


bp N;, pb ONS P N> 

2 17 19 53881 47 9257329 
3 73 23 87481 53 22000801 
5 241 29 117049 59 48473881 
7 1009 31 515761 61 48473881 
11 2641 37 1083289 67 175244281 
13 8089 41 3206641 71 427733329 
17 18001 43 3818929 73 427733329 


79 898716289 
All N, above are primes except for 


Nu = 19-139 
Nu = 47-383 
Ne = 67-1747 


Na = 643 -4987 
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The interest in pseudo-squares is heightened by the fact that they may be 
used in tests for primality, as shown by MARSHALL HALL.‘ The operation of 
the SWAC and the reduction and checking of the output data was done by 
JouHN SELFRIDGE. 

D. H. L. 

1M. Kraitcaik, Recherches sur la Theorie de Nombres, v. 1, Paris, 1924, p. 41-46. 


? D. H. LEuMER, ‘The mechanical combination of linear forms,” Amer. Math. Monthly, 
v. 35, 1928, p. 114-121. 


3D. H. Leumer, “The sieve problem for all-purpose computers,”” MTAC, v. 7, 1953, 
p. 6-14. 


4M. HALL, “Quadratic residues in factorization,’’ Amer. Math. Soc., Bull., v. 39, 1933, 
p. 578-763. 


171.—L. F. RicHarpDsoNn (1881-1953). This English mathematician made 
several notable contributions to numerical analysis. A brief account of his 
life by P. A. SHEPPARD appears in Nature (v. 172, 1953, p. 1127-8). 

His work on numerical analysis (apart from that appearing incidentally 
in his book') was mainly contained in three long papers: 


A: “The approximate arithmetical solution by finite differences of physical 
problems, involving differential equations with an application to the 
stresses in a masonry dam.” Royal. Soc. Phil. Trans., v. 210 A, 1910, p. 
307-357. 

B: “The deferred approach to the limit’—Part I, L. F. RicHarpson, 
Part II, J. A. GAunt, Royal Soc. Phil. Trans., v. 226 A, 1926, p. 299-361. 

C: “A purification method for computing the latent columns of numerical 
matrices and some integrals of differential equations,” Royal Soc. Phil. 
Trans., v. 242 A, 1950, p. 439-491. 


There were minor contributions in 


D: “Theory of the measurement of wind by shooting spheres upward,” 
Royal Soc. Phil. Trans., v. 223 A, 1923, p. 345-361. 


An introduction to some of the material of A, B appeared in 


E: “How to solve differential equations approximately by arithmetic,” 
Math. Gazette v. 12, 1925, p. 415-421. 


and one to some of the material of C in 


F: “A method for computing principal axes,” British Jn. of Psychology, v. 3, 
1950, p. 16-20. 


His work is highly individualistic and his language and symbolism pictur- 


esque. For instance he introduced the terms ‘‘marching’’ problem, for initial 
value problems of the form 


y”’ = ky, (0), (0) given, 
and “‘jury’’ problem for a problem of the form 


ys — Sy® + Sy” —y = dy, y =" =” = 0 for x = + 1. 
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We shall discuss briefly two of Richardson’s contributions to numerical 
analysis. 

Richardson early in his life recognized the importance of the use of 
central difference operators in numerical applications, a fact which had been 
pointed out somewhat earlier by W. F. SHEPPARD. The use of central-differ- 
ence approximation to derivatives suggested that the (local) difference 
between the solution to a discrete problem and that of the continuous prob- 
lem would be a power series in h?, h being the mesh-length. Taking account of 
the first term only, it is possible by solving the discrete problem for two 
values of h, and then eliminating the h?-contribution to obtain a better ap- 
proximate solution. This he called the “deferred approach to the limit.” 
Among his examples were the extrapolation from the perimeter of a square 
and hexagon to that of a circle (see MTAC v. 2 1946, p. 114 and p. 223-4) 
and the evaluation of e. He also indicated the passage to the fundamental 
frequency of a continuous string from that of strings of beads. Richardson is 
mainly concerned with applications of his method to the numerical solution 
of differential equations ; there are, however, discussions in D of a quadrature 
and the solution of a Volterra type integral equation. 

In B he examines in some detail the justification of this process, consider- 
ing two main questions: (1) Are the odd-powers of h always absent? (2) How 
small must / be in order that h?-extrapolation may make an improvement? 
We shall not discuss this paper in detail: Undoubtedly the process is a valu- 
able practical tool, but there are certainly cases where it is unreliable.” 
Richardson was fully aware of the possible failures and difficulties which 
might be encountered in its application and discussed various bad examples. 

The latter part of A is concerned with a detailed study of the stresses in 
a dam, with particular reference to a model of the Assuan dam. 

An interesting remark in A concerns the cost of computation about 1910. 
The unit operation was the evaluation of 


ww + Uw + vs + VE — 400 


and for this the rate was 1/18 pence where m was the number of digits carried. 
An average output in the case » = 3 was of the order of 2,000 correct units 
per week, paying about 28 shillings or about 5 dollars at the then current 
rate. 

The second contribution of Richardson which we shall discuss is an 
algorithm for the solution of a system of linear equations, Ax = b. This is 
to choose an arbitrary x = {x,} and then put 


x(t) = x) 4+ 8 {Ax — b}, 


where the factors 8, are to be chosen. Some suggestions for their choice was 
given in A; an up-to-date study has been given by D. M. Youna.* 

The success of this algorithm can be established by expressing the errors 
in terms of the characteristic values of A. Following Young we assume that 
A has linearly independent characteristic vectors V;,---,V¥, with character- 
istic values \y,---,An. Then the error vector e = x) — x satisfies 


et) — e” ob BAe. 
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n 
If we expand e in terms of the v; as: e& = > c,v; we find 


i=1 


e = Sew: TI (1 + Bd. 


i=1 k=1 


This gives 


Nel = E calla | TL + aaa} < epee 
i=1 k=1 
where 


MO = max | I (1 + pa. 


l<i<n | k=l 


For convergence we have to show that M“ — 0. In practice the sequence of 
B, will often be taken as periodic and then it will! be sufficient if the product 
taken over a period, is less than unity. In this case, if we know, for instance, 
that 

0< <b, r=1,2,---n 


and choose §, so that 0 > 6, >— 2b-', then each factor will be less than 
unity and convergence is assured. 

It is clear that these ideas can also be used in the practical determination 
of the characteristic vectors of a matrix. Suppose we have an approximation 
v to a characteristic vector v; of a matrix A and, for simplicity, assume that 
the only contamination is a component of the characteristic vector Ve. 
Suppose we have v = ¢1V; + C2V2. Let A1,A2 be the characteristic values cor- 
responding to vi,V2. Then, for any a(#),) 


v) = (A wal al)v = oi(1 —_ a)Vi + C2(Aze -~ a)V2 
(Ay —_ a)[o,¥; = Co(Az2 -— a) (Ay sd a)—ve |. 


The strength of the component of vz will therefore be reduced if | (A: — a| 
< |i — al, i.e. if @ is nearer \2 than ),. If this is so successive repetition of 
the multiplication A — al will purify v. This method is a generalization of the 
familiar ‘‘power” method for the determination of the characteristic value of 
largest absolute value. 

In C, Richardson exploits this idea, with a wealth of numerical examples, 
including cases when the matrix is unsymmetric or has non-linear elementary 
divisors. He discusses the use of purifiers of the form 


(A — aI)(A — acl)---(A — al) 


and the optimal choice of the a;. It is clear that information about the 
location of the characteristic roots is essential for satisfactory choice of the 
a; Richardson makes use of the bounds given by Hirscu, Rayleigh’s 
quotient, and the comparison of the ratios of corresponding components of 
Av and v. The latter is used in an intuitive way, no mention being made of 
the result of CoLLatz, that there is always at least one characteristic root 
between the greatest and least of the ratios of the components. 


A 
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1L. F. Ricuarpson, Weather Prediction by Numerical Processes. Cambridge, Eng- 
land, 1923. For the following comment on this, | am indebted to G. E. Forsyrtue, “It is a 
monumental attempt to forecast for six hours, from almost no initial condition, and (I under- 
stand) a poor balance of At and Ax, Ay, Az. It i is superbly written and the author has (in 
my opinion) the most elegant English style of any mathematical writer of the century. 
[See p. 219 of this book, or the first page or two of C.] The Preface speaks for itself of the 
a encountered by the author.” 
f. W. Wasow, “Discrete approximations to elliptic differential equations.” Zeit. f. 
ang. “Math u. Phys., v. 5, 1954. 

*D. M. Youne, “On Richardson’s method for solving linear equations with positive 
definite matrices.” In. Math. and Physics, v. 32, 1953, p. 243-255. Experiments on the 
solution of the Laplace equation by this method, on ORDVAC, have been carried out by 
D. M. Young, and C. H. WaRLIcK. 


CORRIGENDA 


V. 8, p. 93, 1-3, for 124 = p* read 12p + py’. 
V. 8, p. 106, 1. 8, for PEARCY read PEARCEY. 
V. 8, p. 121, 1. 20 for +3(2 + 4) read —3(2 + 4). 


EpiToriAL Norte. With this issue of MTAC the present Editorial Com- 
mittee rounds out its fifth year and resigns. It is a pleasure to thank our 
many contributors, reviewers and referees for their cooperative assistance 
to the Committee and to MTAC. Future editorial correspondence should be 
addressed to 


Dr. C. B. ToMPkKINS 
Numerical Analysis Research 
University of California 

Los Angeles, 24 

California 
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